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Abstract: We study causal wedges associated with a given sub-region in the boundary of 
^ asymptotically AdS spacetimes. Part of our motivation is to better understand the recently 
o I proposed holographic observable, causal holographic information, x, which is given by the area 
fH of a bulk co-dimension two surface lying on the boundary of the causal wedge. It has been 
suggested that x captures the basic amount of information contained in the reduced density 
^ matrix about the bulk geometry. To explore its properties further we examine its behaviour in 
time-dependent situations. As a simple model we focus on null dust collapse in an asymptotically 
00 AdS spacetime, modeled by the Vaidya-AdS geometry. We argue that while x is generically quasi- 
telelogical in time-dependent backgrounds, for suitable choice of sub-regions in conformal field 
O theories, the temporal evolution of x is entirely causal. We comment on the implications of this 
^ observation and more generally on features of causal constructions and contrast our results with 
j> the behaviour of holographic entanglement entropy. Along the way we also derive the rate of 
^ early time growth and late time saturation (to the thermal value) of both x and entanglement 
entropy in these backgrounds. 
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1 Introduction 



One of the important questions in holography is to understand the precise dictionary between the 
bulk spacetime and its avatar in the dual boundary quantum field theory. Over the years we have 
learnt to encode geometry in terms of field theory observables. While there has been considerable 
success in identifying key geometrical features in terms of the field theory data, it is nevertheless 
clear that the translation between the two descriptions is far from complete. We are still trying 
to ascertain the sharpest statement about geometry. The present work, which is exploratory in 
spirit, examines the features of observables generated from purely causal constructs of the bulk 
spacetime. 

One class of questions which probe the CFT encoding of the bulk geometry starts by re- 
stricting the boundary region on which one has access to CFT data. For example, supposing we 
know the full reduced density matrix p_4 on some spatial region A on the boundary, how much 
information does p_4 contain about the bulk geometry? This question was examined recently in 
[1-4], though no consensus on the final answer was reached. Instead of confronting this ques- 
tion head-on, [3] took an indirect approach of asking: what is the most natural (i.e., simplest, 
nontrivial) bulk region associated to Given such a natural and therefore important bulk 
construct, we expect that there should exist a correspondingly important dual quantity in the 
field theory, perhaps waiting to be found. If we succeed in identifying such an object within the 
field theory, we will obtain a more direct handle on the gauge/gravity mapping of the geometry 
and consequently on bulk locality. 

The most immediate geometrical construct associated with A that probably springs into 
the reader's mind is an extremal co-dimension two bulk surface which is anchored at the AdS 
boundary on dA. Indeed, this is a well-known and important construct: In [5, 6], Ryu and 
Takayanagi conjectured that for an equilibrium state, the entanglement entropy of A is given by 
the area of precisely such a bulk surface: a co-dimension two minimal area surface at constant 
time which is homologous to A and anchored on dA (entangling surface). More generally, for 
states that evolve non-trivially in time, one should use extremal surfaces as argued in [7]. Though 
we have no proof to date, there is mounting evidence that entanglement entropy is indeed given 
by the extremal surface area; see [8-12] for arguments to show that the holographic constructions 
satisfy entropy inequalities and [13] for a derivation of the holographic entanglement entropy in 
some special circumstances. 

But is this the most natural bulk construct associated with A7 Finding an extremal surface 
in the bulk requires the knowledge of the bulk geometry. Although this is what we are ultimately 
after, there is a more primal and perhaps more fundamental concept, namely the bulk causal 
structure (the knowledge of which requires a proper subset of the information contained in the 
bulk geometry). In [3] we argued that the simplest and most natural construct is in fact the 
causal wedge (which we will denote 4^) corresponding to the boundary region A, and associated 
quantities. We will review the definition of 4^ in detail below; but to orient the reader, a 
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Fig. 1: A sketch of the causal wedge and associated quantities in planar AdS (left) and global AdS 
(right) in 3 dimensions: in each panel, the region A is represented by the red curve on right, and 
the corresponding surface H_4 by blue curve on left; the causal wedge lies between the AdS 
boundary and the null surfaces 9-|-(#^) (red surface) and (?_(#^) (blue surface). 

succinct description is as follows: Take the boundary domain of dependence O^i of A; this is the 
boundary-spacetime region where the physics is fully determined by the initial conditions at A. 
The bulk causal wedge is the intersection of the causal past and future of (>^- Hence any causal 
curve through the bulk which starts and ends on 0^ must be contained inside the causal wedge 
and conversely we may think of 4^ as consisting of the set of all such curves.^ 

The causal wedge is a (co-dimension zero) spacetime region; but we can immediately identify 
associated lower- dimensional quantities constructed from it, namely bulk co-dimension one null 
surfaces, forming the 'future part' c?+(#^) and 'past part' ^-(♦^) of the boundary of the causal 
wedge, as well as a bulk co-dimension two spacelike surface H_4 lying at their intersection. For 
orientation, these constructs are illustrated in Fig. 1, for planar AdS (left) and global AdS (right). 
Hence, S^, dubbed the causal information surface in [3], is a spacelike surface lying within the 
boundary of the causal wedge which penetrates deepest into the bulk and is anchored on dA. 
In [3, 7] we demonstrated that while S_4 must in fact be a minimal surface within d{^j) that is 
anchored on dA, it is in general not an extremal surface in the full spacetime. There however are 
certain situations where the causal information surface S_4 actually coincides with the extremal 
surface 2^ as noted in [3]. It was conjectured there that the corresponding density matrix p_4 
was maximally entangled with the rest of the field theory degrees of freedom. Below, we will 
consider these special situations further and provide additional evidence for this suggestion. 

So far we have utilized solely the causal structure of the bulk to construct our natural bulk 

^Note that [1] shows that the causal wedge 4^ is equivalently defined in terms of the intersection of future and 
past going hght-sheets emanating from 0^. They further argue using the covariant holographic entropy bounds 
[14] that this implies that the causal wedge 4^ must be the maximal region of the bulk that can be described by 
observables restricted to O^- Since the extremal surfaces computing entanglement entropy necessarily lie outside 
the causal wedge [3, 12] it however seems more natural that the boundary theory restricted to ()a is cognizant of 
a larger part of the bulk as argued in [2] . 
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region 4^ and associated surface S^; but now we finally recourse to geometry. In particular, in 
analogy with entanglement entropy related to the proper area of we identify a quantity 
Xa related to the proper area of E^, 



In [3] we called this quantity the causal holographic information, and studied its properties in 
equilibrium. In particular, we conjectured that xa provides the lower bound on the holographic 
information contained in the boundary region A] however to make this more precise or mean- 
ingful, we need to understand better what sort of quantity x is from the field theory standpoint. 
To that end, we will continue the exploration of the properties and behavior of x under various 
circumstances. 

Since [3] constructed x ^i-nd causal wedge in equilibrium configurations, we will explore 
the properties of x in more general out-of-equilibrium situations. Rather than examining the 
qualitative features of our constructs in arbitrary spacetimes, it will be instructive to obtain more 
detailed quantitative results. To that end, it is useful to pick a specific class of examples, which 
are not only tractable, but also far out of equilibrium. The more extreme the time variation, the 
more easily we can sample the 'dynamics' of 4^, S_4, and xAi e.g. when the spatial position of A 
is fixed but we study it at different times. We focus on a particularly simple time-dependent bulk 
geometry, describing a collapse of a thin null spherical shell to a black hole in AdS, namely the 
Vaidya-AdS spacetime. Since a (sufficiently large) AdS black hole corresponds to a thermal state 
in the field theory, this geometry has been much-used to study thermalization in the field theory 
via black hole formation in the bulk (see [15-25] for a sampling of references). Moreover, since the 
shell is null, the collapse to a black hole (and hence the corresponding boundary thermalization) 
happens maximally quickly. Also, since the shell is thin (and so starts out from the boundary at 
a single instant in time), the change in the boundary corresponding to the introduction of the 
shell is sudden: we deform the boundary Hamiltonian and then let the system equilibrate - in 
other words, such process in an example of a quantum quench.^ We refer the reader to [27-31] 
for further discussions of quantum quenches (including computation of observables) in conformal 
field theories and to [32-43] for discussion of holographic quenches and thermalization. 

To study the behavior of 4^, H^, and xa in the Vaidya-AdS class of geometries, we set 
up the geometrical construction in §2. We discuss the general expectations for the behavior of 
these quantities in §2.1 and derive the explicit equations to construct in §2.2. In subsequent 

^We should nevertheless emphasize that the modeling of quantum quench using the Vaidya-AdS spacetime is 
somewhat contrived. Indeed, in the boundary theory the final state is thermal and known by construction, while 
in a typical global quench protocol one changes a parameter of the Hamiltonian at some time without knowing 
the final state, which is not guaranteed to be thermal. For instance, for simple integrable models like the Ising 
chain it is known that the final state is given by a Generalized Gibbs Ensemble (GGE) [26] where in principle all 
the (infinite) integrals of motion occurs. Furthermore, it is more realistic to introduce localized sources to deform 
the theory, in contrast to the homogeneous disturbance (injected in the UV) used in the null shell collapse. The 
primary advantage of the models we describe below is their tractability. This caveat should be borne in mind 
before drawing general conclusions from our analysis. 
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sections we examine the detailed 'dynamics' of our causal constructs, focusing on the causal 
holographic information xa, for Vaidya-AdSa in §3 and for higher-dimensional Vaidya-AdS in 
§4. While the thin shell Vaidya-AdS class of geometries studied hitherto illustrates many of the 
essential features of the causal wedge and x, some of these results derive from the large amount of 
sjTiimetry which has rendered these examples tractable in the first place. To surmount potential 
bias towards special situations, and in order to gain more intuition on the requirements which 
any putative CFT duals of these quantities must satisfy in general, we comment on some general 
properties of the causal construction in §5. These will be further explicated in a companion paper 
exploring more formal aspects [45]. The appendices collect some of the technical computations for 
three dimensional Vaidya-AdS spacetimes (in particular Appendix B contains some new results 
on entanglement entropy for global collapse and details of early and late time behaviour). 



2 Preliminaries 

As described in §1, we would like to understand the behavior of the causal wedge ♦yi and the 
causal holographic information xa introduced in [3], in situations where the reduced density 
matrix associated with the given spatial region A is time dependent. We will concentrate 
on the process of thermalization following a sudden disturbance which has oft been used as a 
convenient toy model of a quantum quench. In particular, we consider a field theory on globally 
hyperbolic background geometry = Sg x in which we introduce a homogeneous distur- 
bance at an instant in time t = tg. Of specific interest will be the cases where the background 
is either Minkowski spacetime = M'^^^ or the Einstein Static Universe (ESU), T,^ = S'^"^. 
We can generate deformations via explicit (relevant) operators introduced into the Lagrangian, 
and ensure homogeneity by smearing the insertion over the constant time slices Sg. The result- 
ing configuration will then undergo some non-trivial evolution, whose consequences we wish to 
examine for a specified boundary region A C T.^- 

In the gravity dual, the said process of thermalization will be described by a simple spherically 
symmetric null shell collapse geometry. We model this by the Vaidya-AdS^+i spacetime with the 
metric: 

ds' = 2dvdr- f{r, v) dv' + rfSL^,^ , /(r, v) = (^1 + ^ ~ (2.1) 

where r is the bulk radial coordinate such that r = oo corresponds to the boundary, the null co- 
ordinate V coincides with the boundary time (i.e. we fix f = t on the boundary of the spacetime), 
and d'S'^_i describes the metric on a plane (sphere) Sg for = {K = 1), so that K keeps 
track of the spatial curvature of the boundary spacetime geometry. The bulk spacetime (2.1) in- 
terpolates between vacuum AdS^+i and a Schwarzschild-AdS^+i black hole if m{v) — t- {0, mo} for 
V —7- =i=cxD respectively. While any monotonically increasing interpolating function m(f ) will do 
the trick, the simplest examples are obtained in the so-called thin shell limit, when the transition 
is sharp, 

m{v) =moQ{v-ts) , (2.2) 
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with B(x) being the Heaviside step-function. In this case the shell is localized at constant v = tg, 
imploding from the boundary r = oo to the origin r = 0. Moreover, we can write the metric in 
a piecewise-static form, 



where the subscript a stand for i inside the shell and o outside, and the shell separating the 
two spacetime regions is at some radius r = Ra{ta) corresponding to a radial null trajectory. 
Although r is continuous across the shell, the time coordinate t is not. Though these thin shell 
geometries will be main focus of our discussion, we will start by setting up the construction in the 
more general spacetimes (2.1), allowing the deformation to act more smoothly temporally. This 
will enable us to easily derive the jump across the thin shell as well as to check our analytical 
results by numerical computations. 

Having specified the bulk geometry, let us now return to the bulk quantities we wish to 
construct, ^,4, S_4, and xa^ for ^ given boundary region A. It will be useful to start by recalling 
the general story, to better understand the simplification afforded by (2.1) and the choice of 
regions we use below. Following [3], we define^ the causal wedge as 



where the domain of dependence G contains the set of points through which any inex- 
tendible causal boundary curve necessarily intersects A. Both and 4^ are defined as causal 
sets; as such, their boundaries must be null surfaces, generated by null geodesies (within Bd and 
in the d + 1 dimensional bulk spacetime, respectively), except possibly at a set of measure zero 
corresponding to the caustics of these generators. This means that constructing the causal wedge, 
for any boundary region A and in any spacetime, boils down to 'merely' finding null geodesies 
in that spacetime. The crux of the computation typically lies in delineating where these future 
and past null congruences intersect. In practice, though, it is desirable to simplify the prob- 
lem still further, by considering convenient regions A and convenient asymptotically- AdS bulk 
geometries. 

Consider first the construction of 0.A within the boundary spacetime Bd- Although this 
background spacetime is simple (e.g. spherically symmetric around any point), for generic regions 
A, the domain of dependence is as complicated as the shape of A, as it terminates at a set of 
caustic curves, the locus where the generators (namely null geodesies emanating normal to dA) 
intersect. However, for spherically symmetric regions A, the symmetry of the setup guarantees 
that within each (future and past) congruence, all null geodesic generators intersect at a single 
point. 

Hence for any interval in 1 + 1 dimensional boundary or for round ball regions in higher- 
dimensional boundary, the domain of dependence is fully characterized by a pair of boundary 

^ The notation is explained in more detail in [3] . Briefly, by we mean the causal past and future in the full 
bulk geometry, whereas Jg indicates the causal past and future restricted to the boundary. 




(2.3) 



^A = J^[OA]nJ-[OA] 



(2.4) 



- 6 - 



points, corresponding to its future tip and past tip . These two points then hkewise 
characterize the full causal wedge 4^, since definition (2.4) merely extends this construction into 
the bulk, 

0^ = J, [g^] n J+[g^] =^ 4^ = J-[g^] n J+[g^] . (2.5) 

Therefore (9+(#^) C 9J~[g^] and (9_(#^) C dJ^[q^] are generated by null bulk geodesies which 
terminate at or emanate from g^, respectively. 

In general spacetimes, finding null geodesies amounts to solving a set of coupled 2nd order 
nonlinear ODEs, so one would typically need to resort to numerics to construct these. Though 
the equations simplify substantially for spherically symmetric geometries of the form (2.1) as 
presented in §2.2, they still retain the form 2nd order coupled nonlinear ODEs. However, for 
piecewise-static and spherically symmetric geometries (2.3), there are enough constants of motion 
to obtain the geodesies by integration; in fact, in the specific cases of interest, we can even obtain 
analytic expressions for the geodesies. Since the thin shell renders Vaidya-AdS merely piecewise 
static, we need to supplement our expressions for the geodesies in each static piece by a 'refraction' 
law for geodesies passing through the shell. This, however, is easy to derive from the geodesic 
equations for the global geometry, as we explain in §2.2. 

The remainder of this section is organized as follows. In §2.1 we focus primarily on spherical 
regions A and thin shell Vaidya-AdS geometries, to motivate our general expectations for the 
dynamics of the causal wedge and x- §2-2 we go on to derive the equations to calculate these 
constructs explicitly. 

2.1 General expectations for xa 

Consider a spherical entangling region A, specified by its radius a, located at time t = t^. As 
indicated above, one may think of 0^ as enclosed by inverted light cones over the region A, so 
Oyi is equivalently specified by its future and past tip, and g^; clearly for Minkowski or ESU 
boundary geometry, the time at which these tips are located is simply 

tg'^ =t_A-a , tgA =t_A + a . (2.6) 

Armed with this data, we are now ready to describe the qualitative behavior of our bulk constructs 
S^, and Xa thin shell global Vaidya-AdS geometry. 

Let us start by making the following simple observation: If we choose A such that its causal 
wedge ♦yi lies entirely in the AdS part of the geometry, xa will have the same 'vacuum' value as in 
pure AdS. Similarly, if lies entirely in the Schwarzschild-AdS (SAdS) part of the spacetime, 
then Xa will have the 'thermal' value it would have in the corresponding eternal black hole 
geometry. The former will be guaranteed if we take t_4 to precede tg by sufficient amount, such 
that the future tip g^ lies in AdS, namely tgA < tg - then by causality the rest of cannot 
know about the shell. Similarly, for to lie entirely outside the shell, in the SAdS part of the 
spacetime, it suffices that t^v > ts, since then the null rays from g^ can never catch up with the 
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Fig. 2: Radial profile of the causal wedge for fixed = —1.5 (left), = (middle), and = 1.5 
(right), for a set of A, color-coded by size a. The thick black curve on right in each panel is the 
AdS boundary, the dashed black line on left is the origin, the dashed red curve the event horizon 
(whose final size \s = 2 in AdS units), and the thin brown diagonal line the shell. The black 
dots denote the radial position of S_4 corresponding to the given A at time and size a. Our 
coordinates are such that ingoing radial null geodesies are diagonal everywhere (i.e. parallel to the 
shell). The plots are made for Vaidya-AdSs spacetime. 



Regime 


time 


equivalent ly 




04^ a) 


1 


t_4 < — a 


tqV < ts, tgA < ts 


same as in AdS 


same as in AdS 


2 


ts — a < tj\^ < ts 


tq^ < ts, tgA > ts 


same as in AdS 


intersects the shell 


3 


ts < t_A < ts + a 


tq^ < ts, tqA > ts 


intersects the shell 


intersects the shell 


4 


tA> ts + a 


tq^ > ts, tqA > ts 


same as in SAdS 


same as in SAdS 



Table 1: Behavior of boundaries of causal wedge (5_(4^) and d-^{^_A) depending on t_4. (Here 
Schwarzschild-AdS is abbreviated by SAdS.) 

null shell and sample the AdS region. Conversely, if t^v <ts< tq/\, then some part of lies in 
AdS and some part lies in the black hole geometry. To illustrate the point, in Fig. 2 we plot the 
radial profile of the causal wedge for a set of region sizes a for three values of t_4: t_^ < ts — a 
(left), = ts (middle), and tji^> ts + a (right). 

To examine this in bit more detail, in Table 1 we tabulate how the future and past parts 
of the causal wedge boundary behave, depending on tjs,. If both (9_(^^) and d+{^ j\) behave as 
in AdS, then so does S_4 and xa- Similar statement holds for both parts behaving as in SAdS. 
However, the intermediate case has a richer behavior: For < tg, none of the null geodesies 
starting at can cross the shell before being intersected by those ending at q^, which means 
that ^-(♦^) still behaves as it would in AdS. However, despite the fact that Syi lies on this 
surface, it will not be the same curve as in AdS if tgA > tg, since the other null surface 9+ (4^) 
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ends in the SAdS part of the geometry and therefore it no longer behaves as in pure AdS. So 
in the regime tg — a < tjs, < ts indicated in the second hne of Table 1, lies only within the 
AdS part of the geometry, but it is nevertheless deformed from the pure AdS behavior. Since 
the surface S_4 is deformed, one would naturally expect that its area xa is likewise deformed 
from the AdS ('vacuum') value. We will however see later that for the special case of spherical 
entangling regions A., this is in fact not the case; this is one of the surprising revelations of our 
exploration. 

The deformation (from its AdS behavior) of H_4 will grow as t_4 increases from tg — a to tg, 
since more and more of 9+ (4^) samples the SAdS part of the geometry. When t_4 > t^, S_4 itself 
can no longer lie entirely within AdS, since its boundary is in SAdS, by virtue of (92 _4 = OA. 
However, not all of S_4 can be entirely in SAdS either while tgv < tg-, since the radial null geodesic 
from must remain to the past of the shell and hence the deepest part of S_4 remains in AdS 
until = tg + a, when the thermal regime is entered. Since the geometry is continuous, we 
expect Xa to vary continuously (and in fact monotonically) with tj,. 

Hence, the expected behavior of xAi characterized in terms of t^, is: 

< tg — a Xa = vacuum result 

tg — a < tji^ < tg + a Xa has non-trivial temporal variation 

t_A > tg + a Xa = thermal result (2.7) 

This means that by general causality arguments, we expect the following to hold: 

1. The 'thermalization' timescale as characterized by xa scales linearly with the system size.^ 

2. is mildly teleological; it responds in advance to the perturbation on a timescale set by 
light-crossing time of A. 

The fact that xa generically responds to the presence of the shell at an earlier time < tg on 
the boundary follows from the fact our construction involves 0^ which samples both the future 
and past of the boundary region A. While ostensibly peculiar, this teleological nature is capped 
off by the light-crossing time, set by the size of the region. Hence the teleological nature of xa 
is not as bad as it sounds, since if we imagine measuring any thermodynamic quantity which 
pertains to the full system, we would need at least this much time anyway. Of course as the 
system size goes to infinity (in the planar case), we will see the usual teleological behavior often 
associated with black hole horizons. 

Both of these timescales (teleology and thermalization) are given simply by a, which is not 
so surprising since it is the only scale characterizing A. However, we can in fact also generalize 
the above statements to any-shaped region A, with appropriate identification of a: Since the 

^ In other words, here we mean the timescale on which it takes XA to achieve its thermal value after the 
excitation. This is not necessarily the timescale by which all observables in the field theory would achieve their 
thermal values; indeed, depending on the diagnostic we use, we may never see true thermalization. 
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boundary metric is fixed and has a well-defined notion of time t, we can define a as the difference 
between and the earliest time to which reaches (or equivalently, half the timespan of Oyi)- 
Such identification provides a natural notion of characteristic size of the region A, and with this 
definition, statements 1 and 2 above hold for any simply-connected^ region A. In the following 
sections we explore these statements in some detail, starting first with the simplest case of d = 2, 
where we can carry out most of the constructions analytically. 

Before turning to the geodesies which govern our causal constructs, let us make one further 
remark about the nature of x- Above, we have been glibly discussing the 'value' of xa] however, 
the surface stretches out to the boundary of AdS and hence xa is divergent. Moreover, it has 
already been shown in [3] that the divergence structure of the area of S_4 is generically different 
from that encountered in the area of the extremal surface Cyi relevant for the computation of 
the holographic entanglement entropy (though in both cases the leading divergence is given by 
the area law). Hence it is meaningless to compare xa ~ for a given region as a function of 
time, except in special circumstances (e.g. d = 2). We therefore will most often concentrate on 
regulated answer obtained by background subtraction, defining 

SxAit) = XA{t) - ; = lim XA{t) (2.8) 

and similarly for 5Sjs,{t). 



2.2 Geodesies in Vaidya-AdS^+i geometry 

Let us now collect some basic facts about geodesies in the spacetime (2.1) that will prove useful 
in the sequel. Since the full d + 1 dimensional spacetime has spherical (for i^' = 1) or planar (for 
K = {)) symmetry, one can effectively reduce the problem of finding geodesies to 3 dimensional 
problem, characterized by r, f , and (the latter generates a Killing direction of whose norm 
defines our radial coordinate g^^p = r^). Then for an affinely-parameterized geodesic congruences 
with tangent vector = i) d'^ + r d!^ + i^d^, it is convenient to define the 'energy' E, 'angular 
momentum' L, and norm of the tangent vector k: 



E = ~pad^ = fv - r 



K = pap'' = -fv"^ + 2vr + r^ip^ = — - + — + 

f f r 



(2.9) 

^2 -2 ^2 



2 



where ' = Note that we are considering full congruences smeared in the directions orthogonal 
to (9^ in to exploit the symmetry. 

For affinely-parameterized null or spacelike or timelike geodesic, k = or 1 or -1, respectively; 
in particular it is a constant of motion. Since d'^ is a Killing field, L is a conserved along the 
full geodesic. On the other hand, since 9° is not a Killing field, E is in general not conserved. 



We will briefly consider non-simply-connected regions in §5. 
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However, in the thin shell limit it is conserved for each piece of the geodesic (inside the shell and 
outside the shell) individually, which we will exploit. In particular, for constant E, the geodesic 
{f (A) , r(A) , f{X)} can be obtained by integrating (2.9). 

While the first-order equations (2.9) are convenient to use in finding the geodesies analyt- 
ically, we can't solve them globally using integrals when E is not constant. The second-order 
geodesic equations valid for generic f{r,v) are 



V + -f,rV -Tip 







^ + 2 U f,r- f,v) V'^ - f,rrV - r f if"^ 



(2.10) 



(f -\ r (y9 

r 







— ^r,v) etc.. In order to solve (2.10) to obtain a specific 



where we use the shorthand fr — 
geodesic, we need to supply two initial conditions for each of the three coordinates. In terms of 
the initial position {vQ,ro,ipo} and the initial velocity, specified by k, L, and initial energy Eo, 
and also a discrete parameter 77 = ±1 which specifies whether the geodesic is initially ingoing or 
outgoing, these are 



v{0) = vo , viO) 



1 



firo,Vo) 



Eo + V\ + f{ro,vo) 



K 



r(0) = ro 
(f{0) = (fo 



r(0) : 

m 



L2 



K — 



(2.11) 



Usually one can exploit symmetries to set = 0. 
numerically to find any geodesic through the bulk. 



For any given f{r,v), we can solve these 



Though the coordinates {v, r, ip} are useful for finding geodesies, they are not the best for 
visualization since the AdS boundary is at r = cxd and f is a null coordinate; hence on our 
spacetime diagrams (such as right panel of Fig. 1, Fig. 2, as well as many of the following figures) 
we present in this paper, we plot p = arctan r radially and f — p + | vertically, so that ingoing 
radial null geodesies are straight lines at 45°. (Note that except for pure AdS spacetimes, the time 
delay which outgoing radial geodesies experience when climbing out of gravitational potential 
well is manifested by these being generically steeper than 45° lines.) 

Jump across thin shell: We now consider the thin shell limit. Since we can solve (2.9) by 
integration in each part of the spacetime where E is constant, all that remains is to account for 
the jump across the shell. As discussed in [16] (see Appendix F), the jump follows immediately 
from (2.10). With f{r,v) = r"^ + K + Q{v) mo/r'^"^, while / and f^r remain finite with a discrete 
jump, /t, ~ 6{v) diverges ai v = 0.^ Thus for a geodesic crossing the shell, since the coordinates 



^Without loss of generality, we henceforth set tg — 0. 
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{v{\),r{\),(p{\)} are continuous across the shell, {v,r,ip} must also remain finite (as is evident 
from (2.9)). This means that from (2.10), v and (f must remain finite as well, which in turn 
implies that i) and ip are in fact continuous across the shell. On the other hand, f has a 6{v) 
piece from the f^^ term, so r must jump across the shell. We can easily compute this jump by 
direct integration. 



f ~ - /i(r) v'^ S{v) 



r dX 



dv 



1 



/i(r) i) 



V 2 

with /i(r) = nio/r'^''^^ which means that the jump in ^ across the shell is 

dVi dVg 1/y n \ 1 / \ 

:> — = o (/i - /o) = 7) ^^in ■ 

dVi dVo 2 2 
It is however even simpler to read off the jump in E directly from the fact that 



E + r]JE^ + f U- 



1 

is continuous across the shell. A bit of algebra then gives 



L2 



from which we recover 



(/. + fo) E, -vifi- fo) \lEf + f,{K-^] 



Ei — En 



1 



{U-fo)v. 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



We now have all the information required to explore the properties of xa the Vaidya-AdS 
spacetimes explicitly. 

3 Shell collapse in three dimensions 



Having gleaned some general features of xa in time dependent geometries, we now turn to the 
specific example of null shell collapse in AdSs, modeled by (2.1) with d = 2. In this case we take 
dS^ = (iyj^ with the spatial circle parametrized hj ip ip + 27r, and 



f{r,v) 



hi 



for f = t < 
ioi V = t > Q 



(3.1) 



so that the spacetime is global AdSa before the insertion of an operator deformation at = and 
BTZ with horizon radius r/j afterwards. This could be achieved for instance by homogeneously 
injecting energy along the spatial circle. The region A is then taken to be an arc of length 2 y9_4, 
without loss of generality lying between ±y9_4. This region will be taken to lie entirely at constant 
time t = t_4 on the boundary. 

The general strategy for finding the causal wedge is as described in §2. Since the boundary 
spacetime is ESU2 (whose metric is fiat), the domain of dependence of A is given by 0.A = 
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J a b^] '^d i?^] ' ^^"^ correspondingly the causal wedge for A merely extends this construction 
into the bulk, 4^ = J^[g^] fl J^[q^]- Hence to find <9±(#^), and therefore S^, we need to find 
future-directed null geodesies from and past-directed null geodesies in g^. These future and 
past tips lie at 

: = t;^ = t_4 + y9_4 , = , r = cx) , 



(3.2) 



The expressions for the null geodesies themselves in the AdSa part of the geometry are given 
by the following expressions [t] = ±1 for outgoing/ingoing respectively and y^oo = w.l.o.g.): 



f (r) = Voo — — {r] + 1) + 1] arctan ^a/ (1 — k"^) — k"^^ + arctani 

- k^)r^ 



ip{r) = T] 



arctan 



A;2 



k 



sign(fc) 



(3.3) 



with k = L/E for simplicity (using the scaling freedom of the null geodesic affine parameter). 
Likewise, the BTZ null geodesies are (now with j = L/E): 



v{r) =Voo + 



^(r) = r] 



2rh 



In 



_ r^) ( V (1 - f) + f rl-rjrh 



2rh 



In 



a/(1 - f)r^ + frl- jrh 
a/(1 - f)r^+frl + j Th 



(3.4) 



To keep track of various geodesic congruences, it is useful to adopt suggestive^ labels: 



V y-{r^tj and ip /^{r^tj describe outgoing congruence terminating at at the boundary. 
v,^{r,tj and ipn^ {r,(^) describe ingoing congruence starting from . 



In these expressions i stands for the (normalized) angular momentum L/E along the given 
geodesic seqment (which for notational convenience we call k in AdS and j in BTZ); because of 
the refraction (2.15), the value of £ will change between k and j as the geodesic passes through 
the shell. 



3.1 Construction of S_4 

Having the explicit expressions for the geodesies at hand, the desired surface (which is a curve 
m. d = 2) can be found easily. One obvious quantity of interest is the minimal radial position 
attained along the curve S; we denote this by in what follows. It is useful to demarcate our 
discussion into four different time intervals for the temporal location of A, corresponding to the 
four rows of Table 1. We consider these in turn: 

^ Here we envision the boundary as being on the right, as in Fig. 2. 
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1. Vacuum (t^ < —fA)'- Here t/s^,t\j < implying that the entire causal wedge and thus 
is in the AdSs part of the spacetime. Using (3.3) with the initial points (3.2) we chart 

out the surface of i9(#^); see right panel of Fig. 1 for the actual shape when ipj^ = n/S. The 
explicit expressions are unnecessary for our purposes (and can be found in [3]). To determine 
rs, it suffices to consider purely radial geodesies L = k = 0. Equating v />(r, 0) with Vr^{r, 0) we 
immediately find 

r= = cot (y9_4 . (3.5) 
Furthermore, one can conveniently characterize S_4 itself by 

sin p cos ip = cos (3.6) 

where p = arctanr. On our spacetime plots such as Fig. 1, would be a horizontal curve at 
f — p + I = t^. As discussed in [3], in pure AdS (and hence in the present "vacuum" regime of 
Vaidya-AdS), the causal information surface in fact coincides with the extremal surface C^; 
in 3 dimensions this is given by a spacelike geodesic with energy E = and angular momentum 
L = cot ip_A^. In [46] this surface was characterized by 

r\^) = , ^l, . 2 > (3-7) 
cos^ if — sm Lf 

which, as can be easily checked, is equivalent to (3.6). 

2. Shell encounter by 9+ (4,4) only: {—^Pa < tyi < 0): In this time interval, the ingoing 
congruence which generates (9_(#^), specified by^ {f ,^(r, A;_), (^^^^(r, A;_)}, still lies entirely in 
the AdSa geometry as explained in §2.1, cf. Table 1. On the other hand, since > 0, the 
outgoing congruence generating 9+(^^) has segments in both the AdS and the BTZ part of 
the spacetime. Let us denote the segments in the two regions then as {f/,(r, fc+), (y9/,(r, fc+)} 
and {^/.(r, j+), y9/,(r, j+)} respectively, accounting now for the fact that the energies in the two 
spacetimes will differ (while L along an individual geodesic remains constant). 

Starting with the outgoing congruence which terminates at g^, for each outgoing geodesic, 
labeled by j+, we need to find the spacetime point ps = {vs = , , ips} where it hits the shell, 
as well as how does it refract there, specified by the relation between j+ and k^. Using the fact 
that the segment {f^(r, j+), 0/,(r, j+)} connects Ps to we find that 

r, = r, (^coth(r,t;,) + -^,==csch(r,t;,)j , e '^^^ = ^^^^^ ^^^^-^ ^ ^^^-^ . (3.8) 

With the knowledge of r^, we can then solve the refraction condition (2.15) with = L/Eo and 
k^ = L / Ei to find that 

f. ^ 2i+r, (rg-r^) 

^ (2 r2 + 1 - rl) + {rl + 1)^(1" jI) rl + jI rf ' 



^ We now distinguish the angular momenta characterizing the top and bottom of the causal wedge d±(^jC) by 
subscript k± for AdS and j± in BTZ. 
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where Vg itself depends on j+ as given by (3.8). 

The main distinguishing feature of the time interval under focus is that A;_|_(j+) from (3.9) 
spans the entire range ±1. This in turn implies that we can view as the data characterizing 
the full angular span of (9+(#^), and confirms that S_4 lies entirely in the AdSa part of the 
spacetime. 

Having described how the outgoing congruence refracts through the shell, it only remains 
to find where it intersects with the ingoing congruence emanating from g^. For each pair of 
intersecting geodesies, we denote their intersection by = {v^ , , ip^}, which we can determine 
by solving 



While the expressions themselves are easy to write down and solve explicitly as we describe in 
Appendix A, it is convenient to solve (3.10) numerically to find S_4. Note that (3.10) gives a one- 
parameter family of solutions for (with corresponding angular momenta k±), which determines 
S_4. We can naturally take H_4 to be parameterized by j+ G (—1,1), or more conveniently by 
f^x £ {i^E,oo) (for each half of S). Since /c+ = when j+ = 0, we can easily find the minimal 
radial position rs attained by analytically: 



Note that in the relevant regime, rs is a monotonically increasing function of (for fixed v?^ 
and r/j). 

In the left panel of Fig. 3 we plot S_4 (thick blue curve) for a region A (thick red curve), along 
with representative generators of d±{4A) (thin null curves, color-coded by r^), for (y9_4 = 2tt/5 
and the final black hole size rh = 2. (Hence the radial null geodesies drawn as thin red curves are 
precisely analogous to the curves demarcating the causal wedge profile in Fig. 2.) For comparison, 
we also show the extremal surface (thick purple curve). We see that unlike the previous case, 
in this regime is no longer plotted as purely horizontal curve, but rather bends outward and 
downward - i.e. to the past of the constant t = f — p + 1 surface. On the other hand, the extremal 
surface 2yi remains undeformed since, being anchored at t_A < tg = 0, it cannot yet 'know' about 
the shell. While the downward bend of H_4 is easy to see in Fig. 3, the outward deformation is 
more apparent from when viewed from a different angle. To that end, in Fig. 4 we present the 
same constructs as in Fig. 3, but viewed from top, i.e. projected onto a constant time slice. This 
projection is known as Poincare disk, where p is the radial coordinate and ip the angular one. 
Here it is evident that lies closer to the boundary than For orientation we also show 
the final black hole size r^, even though the generators are not directly dependent on it. On the 
other hand, we can see that the generators of d^{^X) are refracted by the shell (which we don't 
show since its projection covers the full Poincare disk and each generator intersects it different 
time and radial position). 



t;K^(rx,/c-) = V: 



ipy,{r^, k+) = (pr^{r^, k-) = if: 



(3.10) 




(3.11) 



15 




Fig. 3: A plot of the causal information surface H_4 (thick blue curve) along with representative generators 
of d±{^_A) (thin null curves, color-coded by r^), in Regime 2 (left) and 3 (right), as discussed in 
text. For orientation we also show the boundary, imploding shell, corresponding event horizon whose 
final size is rh = 2, the region A (thick red curve) whose size is = ^ and time t_4 = —0.1 
(left) and t_A = 0.6 (right), the corresponding domain of dependence 0^ (thin grey curves) with 
its future and past tips q^, as marked, as well as the extremal surface C;^ (thick purple curve) 
for comparison. 




Fig. 4: Top view of the same plot as in Fig. 3 (with the same color-coding scheme), i.e. all curves are 
projected onto the Poincare disk. For orientation, we also indicate the final black hole size 
(dashed red curve). 

3. Shell encounter by S_4: (0 < < fjC)' We now come to the most complicated regime 
of interest (cf. 3rd line of Table 1). For the time interval under consideration, while we still 
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have f V < and f a > 0, there is a quahtative change in the behavior. This is because 
itself will intersect the shell, which one can argue for as follows. Along the ingoing congruence 
{vK^{r,k_),(f)^{r,k_)}, the radial null geodesic (A;„ = 0) stays at f = f v and thus is parallel to 
the shell's trajectory and has to remain in AdS. On the other hand, the geodesies with k_ ^ ±1 
lie close to the boundary and these must intersect the congruence from at f = > on 
the boundary. The only way for this to happen is for the ingoing congruence itself to cross over 
through the shell and sample regions with both signs of v. 

Hence in this regime, both the congruences generating (9+(4_4) and (9_(4_4) refract through 
the shell. The transition (determined by where S_4 itself intersects the shell) is given by some 
critical angular momenta and fcl demarcating the transfer of refraction from the top boundary 
of the causal wedge to its bottom boundary. We have to determine these to find S^. 

The analysis however is straightforward; start with the radial geodesies for which only the 
refraction of the j+ = k_^. = geodesic matters. This is of course similar to what we encountered 
in Regime 2 and it hence follows that the minimal radial position attained along S_4 continues 
to be given by (3.11). We then increase j+ and follow S along its path through the AdS region 
as before. At the same time we monitor the ingoing geodesies along d-{^X) and ask when they 
hit the shell. This happens for 



1 /7 \ / cosfv — Wl — ki . 

Ts = cotvy , cscfv , V's = TT sign(fc_) + arctan . (3.12) 

a/1 - F 2 y k- smvy 

The critical angular momenta j*^ and /c* at which crosses the shell is then obtained by equating 
Ts and !fs in (3.12) with the corresponding result in the black hole part (3.8). Denoting the 
spacetime point where these critical geodesies with j"^ and k^ intersect (which is simultaneously 
the point where S_4 intersects the shell) by px = {vx , fx , Vx}, we have vx = 0, rx = rs{j\^) = 
r,{k*_) and (^x = '^s{3*+) = ^s{kl). 

The strategy for finding then is similar to what was employed in (3.10). For Ij+I < 
or equivalently for r= < r < r^, the previous analysis carries over unchanged. For larger values 
of angular momenta (r > r^), we must first account for the refraction of the ingoing congruence 
from through the shell, by employing the relation between fc_ and j_, analogous to (3.9) and 
obtained from the same refraction condition (2.15), now with rj = —1, j_ = L/Eq and k^ = L/Ef. 

2 k_ r, (r^ + 1) ^ ^ 

^ ' ^ (3.13) 



r, (2r2 + l-rl) + {rl + 1) ^{l - A:^) - kl ' 

The analog of (3.10) which we need is simply obtained by replacing A;_ — )■ j_ and A;+ — )■ j+ since 
the intersection happens in the BTZ part of the spacetime now. 

In the right panel of Fig. 3 we plot along with representative generators of i9-|-(^_4) for 
this regime, as well as the extremal surface €^ for comparison. We can see that S_4 now deforms 
to an even larger extent than in Regime 2 (cf. the left panel), being pushed further outward and 
downward, as well as kinked by the shell. The behavior of the extremal surface is likewise 
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more complicated than in the previous two cases (whether or not crosses the shell depends on 
the interplay of and ^9^; in the present case it does), and its detailed structure will be presented 
elsewhere [47]. However, we can make the general statement that does not coincide with S_4, 
and reaches deeper into the bulk, as characterized by the minimum radius attained along (£_4, 
rg < th- We will revisit this point in the Discussion. 

4. Thermal (t_4 > (fjC)- Finally, let us consider the regime V/\,V\j > 0, so the entire causal 
wedge is in the BTZ part of the geometry. As demonstrated in [3] , the causal information surface 
S_4 now again coincides with the extremal surface (E^, both are deformed outward and downward 
by the presence of the black hole, such that = > r^- By similar arguments as for regime 1, 
we now find the minimal radius reached to be^ 

rs = Th coth (r/j ipj^) = . (3.14) 

The static case expressions (3.5) and (3.14) are in fact the special limits of (3.11) as t_4 — ?■ ±(/9_4, 
respectively. As remarked above, r= increases monotonically with t_4, so in particular the thermal 
result (3.14) is larger than the vacuum result (3.5). 

Now that we have covered all 4 qualitatively distinct regimes, we summarize our results. In 
the left panel of Fig. 5 we plot the actual surfaces S_4 (now color-coded by tX) as varies across 
the 4 regimes, again for a fixed value of = 27r/5 and the final black hole size = 2 (so the 
thick blue curves in Fig. 3 are specific examples of these). We present the same curves S_4 
both on a spacetime plot (left) as well as its projection onto the Poincare disk (right). We can 
clearly see how the surfaces deform outward and downward so as to remain outside the event 
horizon. The 4 regimes are demarcated by the regions A for t_4 = —fA^^jfA labeled in the 
left panel, and we can see that in regimes 1 and 4 the shape of S_4 remains the same, while in 
regimes 2 and 3 the shape of S_4 changes with as expected. The qualitative difference between 
the latter two regimes can be seen if we shift all S's such that they are anchored at the same 
position on the boundary. Then one can confirm that in regime 2, all S's lie on the same null 
surface, while in regime 3 they don't. 

To characterize the change in S_4 under variations of t^, it is better to concentrate on one 
salient feature of S_4 rather than its entire shape. One such handy quantity is the bulk depth 
to which H_4 penetrates. In Fig. 6 we plot the minimum radius (blue curve) and rg (purple 
curve) attained by the causal information surface and the extremal surface respectively, 
as varies across the 4 different regimes discussed above, again for = 27r/5 and final black 
hole size Vh = 2. We clearly see that the expectations explained in §2.1 pan out: rs coincides 
with in regimes 1 (AdS) and 4 (BTZ) and differs in regimes 2 and 3 (shell encounter); in 
particular > T|J (i.e. doesn't penetrate as deep into the bulk) in the latter cases. Moreover, in 
regime 2 (t_4 < tg = 0), while r,£ remains at its AdS value by causality, already starts to vary, 
illustrating the quasi-teleological behavior of 4^. 

^ Note that we have denote the minimal radius attained by and £^ in BTZ spacetimes as for future 
convenience. 
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Fig. 5: (Left): The curves H_4 (color-coded by i^) for a range of t_4 sampling across the 4 regimes 
(separated by the three transitions at = — (/9^,0,(/?^ as labeled; the thin gray curves represent 
A at those times) in increments of O.lc/?^, for ip_A = 2tt/5 and r/j = 2. (Right): Same curves H_4 
projected onto the Poincare disk, analogous to Fig. 4. 




Fig. 6: Comparison of minimum radii (blue curve) and r(C (purple curve) attained by the causal 
information surface H_4 and the extremal surface C;^, respectively, as a function of t^, for the same 
parameters as in Fig. 5, namely (/?_4 = 2tt/5 and r/j = 2. The regimes 1,2,3,4 are again demarcated 

by tA = -tpA,0,ipA- 



3.2 Determining xa smd comparison with 

Now that we have analysed how the causal wedge and the corresponding causal informa- 
tion surface 'evolves' during a thin shell collapse, let us turn to its proper area, the causal 
holographic information xa- particular, we would like to compare the regulated value of xa 
with the regulated entanglement entropy S^- One might naively expect that xa and S^. would 
evolve with in a manner which is qualitatively analogous to that of and rg plotted in Fig. 6. 
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cpA = 1.256 cpA = 0.85 

Fig. 7: The variation of XA with time tj(. We plot the absolute value of XA evaluated with a radial cut-off 
Too = 100 in the Vaidya-AdSa spacetime. The shell implodes from the boundary at tg = 0. The 
plots are shown for two different region sizes indicated above and for different final black hole size 
for each choice of (^_4. 



Indeed, we find that in regimes 1 and 4 (AdS and BTZ), xa and 5*^ must coincide, since the 
actual surfaces whose area these quantities measure also coincide. 

In particular, as can be verified by explicit computation, in the AdS ('vacuum') case,^° 

CcfT 



Regime 1 : 
while in the BTZ ('thermal') case. 



Xa = Sa= ^ log [2 Too sinipj] = ^Ads(v5^) 



(3.15) 



Regime 4 : 



Xa = S. 



Ceff 



A 



log 



2r« 
rh 



sinh {rh (Pa) 



= ^BTziv A, rh) 



(3.16) 



where Too is the radial UV cut-off to regulate the standard divergence encountered in the expres- 
sions. We have also introduced new definitions for the values of the x S in the AdS and 
BTZ geometries respectively for future convenience. 

In the intermediate regime (Regimes 2 and 3) where the causal wedge encounters the shell 
and no longer coincides with we can compute xa numerically. (If fact, in the present 
case we can also use a trick, explained in §3.2.1, to obtain xa almost analytically.) Using the 
more obvious numerical method, we integrate the length element induced from (2.1) onto the 
curve S. This boils down to integrating the proper length in AdS for < r < rx (which for 
Regime 2 is rx = C)0 so this gives the full answer) and using the BTZ metric for r > rx- We 
regulate the result by integrating the length element up to r = Too; since at the end of the day 
we are going to use background subtraction, all we need to do is to ensure that we pick the same 
UV regulator for the AdS spacetime. 

In Fig. 7 we plot the behavior of xa two different values of (pj^. While we see that xa 

^°The expressions are written in terms of the field theory central charge Ccs which is related to the gravitational 



Newton's constant via the standard Brown-Henneaux result 



Ccff 



3 RAdS 



2G 



(3) 
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indeed monotonically interpolates between the AdS value and the BTZ value, we encounter a 
surprise: xa behaves causally: it does not start to grow till the shell encounter! In other words, 
Xa = *S'_4 even in Regime 2, despite the fact that the surfaces H_4 and differ. 

The reason that xa remains at the AdS value in Regime 2, and only starts to vary in Regime 
3, is the following. As explained in §2.1 (cf. Table 1), ^-(♦^) lies entirely in AdS, so lies on 
the same null surface dJ~^[q'^] as the extremal surface 2^ (the latter lying on dJ~^[q'^] by virtue of 
coinciding with the causal holographic surface in pure AdS). In particular, past-directed outgoing 
null geodesies emanating in a normal direction to CB_4 thus generate c?J^[g^], with S_4 lying along 
these generators between (£_4 and g^. 

Since 9J+[g^] is a boundary of a causal set, its generators must be null geodesies which reach 
the boundary at without encountering caustics along the way. By Raychaudhuri's equation, 
this in turn implies that these generators cannot contract towards the boundary, i.e. that their 
expansion along past-directed (outgoing) direction must be non-negative, but cannot increase. 
On the other hand, as shown in [7], the extremal surface is precisely the one with null normal 
congruences (both ingoing and outgoing ones, or both future and past-directed ones) having zero 
expansion. Since the generators of dJ~^[q'^] start out at with zero expansion towards the 
boundary, they have to maintain zero expansion all along the entire 9J^[g^], which can be also 
checked by explicit calculation (as we do in §4). 

Having established zero expansion along null generators of (9J^[g^], the final ingredient in 
our argument is translating this into comparison of areas of and S_4. Since the expansion 
is the differential increase in area along the 'wavefront' of these generators, = ^logy4(A), if 
0(A) = 0, then the 'wavefront' area A{X) stays constant. Furthermore, since we can think of 
as lying at A = and S_4 as lying at A = Ah > Ag (using, if necessary, the freedom of overall 
rescaling of affine parameter along a null geodesic, and noting that at the boundary, finite A flow 
degenerates to a point, so that all constant A wavefronts remain pinned to dA), the increase of 
area between (£_4 and S_4 must vanish, i.e. xa = Sa- 

Note that the above argument crucially relied on the fact that G 97+ [g^]. This situation 
is general for d = 2 where our region A is just an interval, since then 0^ is specified by the 
two points g^ and g^ for any A. On the other hand, as pointed out in §2, this clearly does 
not hold in higher dimensions for generic shapes of A. It is only for special (round) regions 
A that coincides with S_4 in AdS and hence can lie on 9J+[g^]. For generic (non-round) 
A, (Ba does not lie on 9J+[g^], so we do not have a handy curve on 9J+[g^] on which we are 
guaranteed to have zero expansion. On the other hand, this lack of proof does not necessarily 
imply inequality between Sa and xa a-priori. To see whether xa does behave teleologically as 
expected, or whether it still maintains causality for a more subtle reason, in §4 we examine these 
quantities explicitly in higher-dimensional thin-shell Vaidya-AdS for both round and non-round 
regions. 
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3.2.1 Trick to evaluate x^(t^) 



Above we have argued that in the thin shell Vaidya-AdSs spacetime, xa behaves causally, in 
the sense that it stays at the AdS value for all < 0, i.e. up till the appearance of the shell. 
However, it is also clear that between t_4 = and = (when xa saturates to the BTZ value 
=^BTz), there is a non-trivial variation in xi^A) which we evaluated numerically; cf. Fig. 7. The 
numerical computation follows the logic outlined earlier to find S_4 and then evaluating its area 
(further details are presented in Appendix A). 

In the present special case of thin shell Vaidya-AdSa we however can exploit a trick to give 
a simple compact expression for xa which only uses the critical angular momenta j"^ and /c* 
discussed above (3.12) and the forms of null and (zero-energy) spacelike geodesies in pure AdS 
and BTZ. While the expressions for and /c* require solving transcendental equations, we know 
analytically the expressions for null and spacelike geodesies in AdSs and BTZ spacetimes, which 
suffices to bring the expression for xa iiito a convenient compact form. 

The basic idea is simply an extension of the one used to argue that xi^A < 0) = x(tA < —Va), 
now applied to light cones in both AdS and BTZ. Fix t_4 G (0, <fX), and consider the curve = S 
(we drop the subscript A for the time being). This is composed of a central piece which resides 
in AdS and the edge pieces which reside in BTZ and these intersect at Px = {v = 0,rx, ±ipx}- 
In fact, since everything is reflection-symmetric around = 0, it will be convenient to deal with 
only one side (say for positive ip); we'll denote the respective parts of one half of the curve by HAds 
and Sbtz respectively, and correspondingly their proper lengths by ^Ads and Cbtz, respectively. 
The total length of S then determines oc 2 (Z^Ads + -Cbtz)- 

Since the segments ^Ads and Sbtz lie in different geometries, it is useful to tackle them 
separately. A-priori to compute the two contributions to x we would be satisfied with any 
mechanism for computing the respective lengths without actually knowing the form of the curves 
themselves. While this is usually tricky, in the present case we can map the computation of the 
lengths Z^Ads and C-qtz, to computations to lengths of two other known curves in the AdS and 
BTZ spacetime. 

Imagine cutting (half of) S at the intersection with the shell into its two segments at px- 
Since S lies on the boundary of the causal wedge, it follows that ^Ads hes on (9J+[g^]; similarly 
Sbtz lies on dJ~ [q^] . We are going to slide the two segments along these light cones to a point 
where we encounter some known curves whose length is easy to compute. 

First we however have to understand why we are free to slide the curve along the light cones. 
Consider the AdS part of H: by construction, ^Ads lies in the AdS part of the geometry, on 
the light cone (9_(^^) G 9J+[g^], whose null generators have zero expansion (as argued above). 
This means that the length of 3 a is the same as the length of any other curve on dJ^ [g^] which 
traverses the same set of generators, namely those null generators of (9J^[g^] with sub-critical 
angular momentum k < k*_ (we take by convention fcl > on E^Ads without loss of generality). 
So as long as we slide SAds up by the same affine time along the generators of dJ^ [q^] with the 
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ends on the generator A;l we won't change its length. 

Among curves that he on this hght cone in AdS, a particularly convenient one is the zero- 
energy spacelike geodesic in pure AdS, S_4. We know it coincides with S in AdS and therefore lies 
on the future light cone from g^, and moreover its length is given simply by its affine parameter 
which is easy to evaluate. Such spacelike geodesic would lie at constant t in AdS, and in 
particular encounters the critical null generator (i.e. one from with angular momentum fc* ) 
at t = tj[. So £Ads is given by the affine parameter A(r) of a spacelike E = geodesic in AdS 
anchored on dA (i.e. with angular momentum L = = cot v^^), evaluated at the value of r at 
which this geodesic intersects the critical null generator with k = k*_. 

Using explicit expressions for the null geodesies in AdS (3.3), we learn that the element of 
the ingoing congruence with angular momentum k*__ starting from makes it to t = tj, at a 
radial position = l_Q,r^2 ■ Then it is a simple matter to integrate spacelike geodesic r(A) to 

infer \{r^). We use ^ = ^ (r^ + 1) (r^ — r^) for a -E = spacelike geodesic with L replaced by 
the minimal radial position attained. The integral is simple to evaluate and we obtain 

^Ads = 2 log . (3.17) 

Note that (3.17) is independent of y^^, which it has to be by scaling invariance of AdS. Also 
note that for tyfi < 0, we have /cl = 1, so when the entire S lies in AdS, we recover the usual 
divergence in its length. 

Let us now turn to the outer piece of S, namely Hbtz- This part of S lies entirely in the BTZ 
part of the geometry, and in particular on the past light cone 9+(#^) G 9J~[g^]. We again slide 
this down to a convenient position staying on this light cone; the main difference is that we are 
interested only in the segment of the light cone generated by null geodesies with — 1 < j+ < 
with indicating being the anchor point of our slide (having chosen positive k_ we now need 
to choose negative j+). 

Since in pure BTZ, an extremal surface €^ also coincides with the causal information surface 
S_4, the generators of dJ~[q^] must have zero expansion everywhere in BTZ, by the same type 
of argument as for the AdS light cones: S forces the generators to start with zero expansion, 
and Raychaudhuri equation ensures that the subsequent expansion does not grow and does not 
become negative - i.e. it has to stay zero. So it then follows that the length £btz of Sbtz is the 
same as the length of any other curve on dJ~[q^] which traverses the same generators, in this 
case characterized by super-critical angular momentum, j+ < j*^. 

The calculation then proceeds exactly as above; we can pick the spacelike £" = geodesic 
in pure BTZ geometry ending at t = t_4, and find where it intersects with the null generator of 
the past light cone from q^ with angular momentum j*^. Using (3.4) for the explicit form of the 
null geodesies in BTZ we find that rl = (r| — (j^)^)/(l — (j+)^).^^ Integrating the expression 

^^Note that since we are moving the segments of the curves SacIS and Sbtz the radial positions r.^, in AdS and 
BTZ are unrelated. 
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for the spacelike geodesic with E = and L = (again set by the minimal radius attained) 



which takes the simple form ^ = yi^"^ — {'f'^ — f'^) in the BTZ spacetime, between r* and 
the radial cut-off Too, we learn that 



Prom these two simple expressions it follows that in the regime < < 9?^ we have a 
compact expression for xa 



where we have written the expression in terms of the BTZ value of x cf-, (3.16). Note that fc* > 
and j!^ < additive logarithmic piece can a-priori be positive or negative. However, since the 
presence of the black hole effectively repels geodesies, > \k^\, which forces the second term 
in (3.19) to be negative. Moreover, explicit numerical solutions for the critical angular momenta 
confirm that is always negative and xa < -5^btz which is consistent with the numerical results 
of Fig. 7. We would like to emphasize that this is a-priori rather remarkable since the surface 
S_4 lies nowhere near any extremal surface in the bulk, as is evident from Fig. 3. Despite the 
apparent non-locality in the definition of the causal holographic information, the final result is 
manifestly local. We will return to this point in §5. 

3.2.2 The behaviour of xa - Sa 

Having understood how to compute xa, let us finally consider the difference between xa and 
in regime 3, which is the only domain in t where it is different from zero. 

First of all we recall that xa a^nd have a leading area law divergence in the UV which 
is replaced by the logarithmic behaviour in d = 2. Unlike the higher dimensional examples, 
here neither has any further divergences, so it makes sense to consider the difference xa ~ 
in the present case. One naively expects [3] that in this regime xa > S^, since the surface S_4 
lies closer to the boundary and hence ought to have greater (unregulated) length. However, 
it is clear that this cannot be the entire story since we have argued that xa = in regime 
4. It therefore must follow that the difference xa ~^ is non-monotonic. Indeed, our explicit 
computation bears this expectation out. In Fig. 8 we plot variation of xa ~ with time t_A- 
We see that this vanishes at both endpoints of this regime, = and t^, = v^^, and is positive 
in between. Moreover, the slope vanishes at both ends (though the numerics are not well under 
control there). To get a handle on the behaviour of xa ~ near = and = fA we turn 
to an examination of the two quantities in these regimes in a perturbation expansion in time. 

The behaviour for — 0: Firstly, consider the behaviour near = 0. As we explain in 
Appendix B it is quite straightforward to work out the rate of growth of 5*^ from the vacuum 

^^The divergent logarithmic contribution conies from the fact that the curves approach the boundary normahy. 
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Fig. 8: The variation of XA ~ Sa with time for the Vaidya-AdSs spacetime. The plots are presented 
for different values of region size ip^ and final black hole radius for comparison. These have 
been obtained by using (3.19) and (B.31) which are evaulated numerically. We note that the result 
is in good agreement with data obtained by numerically solving for H_4 as described in Appendix 
A and thence computing XA and similary for 
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indicating a quadratic growth in the holographic entanglement entropy about its vacuum value. 

The behaviour of xa can be computed directly using (3.19) which is a clean local formula. 
Were it not for this it would be quite hard to estimate the change in xa about its vacuum value. 
We solve (3.8) and (3.12) for and k*_ for small t_4, which can be done analytically; plugging 
the result into (3.19) we have (with rg = coty9^ and defined in (3.14)): 



XA{tA, ^a) = ^AdS + 
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From (3.20) and (3.21) we conclude that the leading and first subleading terms in the growth of 
Xa and 5*^ cancel each other off leaving behind a cubic growth: 



(3.21) 



Xa 



S 



A 



Ccff 

24 



(cot - Th cothr/iV5^)(r^ csch^ r^LpA - csc^ Lfj, 



r 



h 



l)t\ + 



(3.22) 



^■^As far as we are aware this is the first analytic result in the literature regarding the rate of growth of Sa. at 
early times for finite region size. The linear behaviour in the intermediate times has been noted before since this 
matches the CFT computation quite nicely. We also note that earlier [44] derived an universal expression for the 
early-time growth focussing on arbitrarily large regions in the context of field theories in M^^^. More specifically 
our results are valid for Ia ^ {^Ht'^a}, with no hierarchy implied between the thermal scale set by and the 
region size i^j^, whereas consideration of arbitrarily large regions requires <C r/j ^ (pA- The latter is only 
sensitive to the IR part of the entanglement entropy and does not for example see the saturation to the late time 
thermal value as we describe next. We thank Esperanza Lopez for discussions on this issue. For completeness we 
present in Appendix B the general behaviour of the growth of S'^(i) at early times starting from various initial 
configurations (global or Poincare vacuum and thermal state) cf., (B.32) and (B.34). 
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The coefficient of the cubic is positive definite guaranteeing that xa > in the neighbourhood 
of the origin. 



The behaviour for t_4 ~ (/?_4: At the other end where approaches y?^, the quantities xa 
and Sj[ tend to their BTZ values ^btz respectively. We can use a perturbation expansion in 
e = — to figure out the rate of approach. For Sj^ this is described in Appendix B; the 
upshot of the computation is that approaches the thermal value as a power law with leading 
exponent |. Specifically, 



SA{tA,'PA) = ^BTZ 
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Likewise we can use the result (3.19) to figure out the behaviour of xa in this regime. The 
strategy involves solving (3.8) and (3.12) for j*^ and perturbatively in (y9_4 — t_4 and then 
plugging the result back into the expression for A straightforward algebraic exercise leads to 
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From (3.24) and (3.23) it follows that 
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6 3r2coth^(r;,v5^) 
implying that the curves in Fig. 8 approach the axis quadratically. 

We should note that in the vicinity of = and = (y9_4 the difference between xa and 
S'_4 is smaller than would be anticipated. In the former regime whilst both deviate from their 
vacuum value quadratically, the leading deviations cancel and the cubic term in xa dominates. 
On the other hand for t_4 — )■ ip^ it is the quadratic term in 5*^ which gives the rate of approach 
to the thermal answer with the leading | power canceling out. We would like to suggest that the 
smallness of the difference between xa and S_a has to do with the specific nature of entanglement 
in 1 + 1 dimensional CFTs, a point we will return to in §5. 

Let us also note that from the numerical analysis we see that the difference xa ~ Sj[ has a 
characteristic peak, which lies in the vicinity of ~ | y^^. The location of the peak is mildly 
dependent on both the black hole size and the size of the interval. 



4 Shell collapse in higher dimensions 

Let us now turn to the higher dimensional examples. Here we have a richer set of choices for 
the shape of the region A. We will however restrict attention to two simple examples (disks 
and strips) to illustrate the basic features of the causal wedges and xa- The choice of regions is 
dictated both by tractability and to motivate the general lessons about the causal construction 
one can infer from them. 
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(i). Spherical entangling region: We choose either a ball shaped round region A C M"'^^ 
in Poincare AdS, or slices of the boundary sphere at constant latitude A C S'^^^ in global AdS, 
depending on whether we want to consider field theories on Minkowski space or on the Einstein 
Static Universe. 

For such spherical entangling regions it was argued in [3] that the causal information surface 
coincides with the extremal (in fact minimal) surface S^. So in the vacuum we expect that 
Xa = Sa- However, this ceases to be true once we excite the state. In the process of the 
collapse, which we continue to model by the Vaidya-AdS^+i geometry, we expect to see both 
Xa ^iid S'_4 increase monotonically from their vacuum values. For S'_4, which evolves causally 
with 6Sa = for tA ^ this was seen originally in [7] and has been throughly explored in the 
recent investigations of holographic quench scenarios mentioned in the Introduction. Explicit 
computations confirm a similar result to hold for xa- particular, in the regime < the 
difference 6xa = implying the xa also behaves causally for the same reason as for the d = 3 
case described in the previous section. 

One can in fact prove this analytically; for completeness and to bolster our arguments in 
§3.2, we take a brief detour to show why the area of S_4 remains unchanged for all tA < 0. As 
explained above, in this regime, the surface S_4 lies entirely in AdS, and moreover lies on the 
same light cone dJ~^[q'^] as Without loss of generality, consider AdS^+i in static coordinates 

ds^ = 5 '^-^ 4.1 

where f is the boundary radial coordinate, and z is the bulk Fefferman-Graham radial coordinate. 
Let the d — 1-dimensional region A on the boundary z = be at t = a, f = a, so that the light 
cone in question (corresponding to <9J"*"[g^]) is simply the one from origin, described by 

- + f2 + = . (4.2) 

As a co-dimension 1 surface, the light cone is parameterized by any two of these three coordinates 
and all the angles in fid_2, which just come along for the ride. Now, since both the metric as 
well as A is spherically symmetric in the so will be any putative surface S_4; such spacelike 

d — 1 dimensional surface will then be parameterized by e.g. z (or f or t) and the d — 2 angles, 
but it won't depend on the angles. That means that is specified by a single function, f{z). 
It can be drawn as a curve in {t, f, z} space, with t{z) determined from (4.2). We will now show 
that the area of any surface \& specified by arbitrary function f{z) is in fact independent of f(z), 
and indeed even independent of f(0) = a. The area of ^ given by 



A^ = Vn dz V-t'{zy + f'{zf + l = Vn dz — ^7=^ , (4.3) 

Jo z'^ ^ Jo V + z^ 

where Vq = 27r^/r(^^) is the area of unit S^^"^, z^i^^ is the maximal reach of \E' and the 
last term was obtained by using (4.2) to simpUfy t'{z). We can now use the change of variables 
u = — , to rewrite (4.3) as 

foo d-2 

A^ = Vn , ^ du. (4.4) 

Jo VM + 1 
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Note that not only the integrand, but also the limits of integration, are independent of the shape 
of the lower limit follows from f(zmax) = and the upper limit from z = while f(0) = a > 0. 
The expression (4.4) is useful to obtain the actual area in d + 1 dimensions. 

Note that the total area Ai^, is of course infinite, so we regulate the expression using a finite 
z cutoff. Then slight care is needed to correctly compare the regulated areas of different surfaces 

A consistent cutoff must be one which is mapped between different surfaces by the null 
generators of the light cone, and is implemented by keeping Uoo = t{zco)/ z^o fixed. For example, 
if we fix the cutoff e for the original surface given by f(f{z) = y/a"^ — z"^, then along another 
surface \E' specified by f^(;z) ~ a + f'^(O) 2;^ + . . ., the new cutoff is modified at quadratic order. 



However for t_4 > we expect to find xa 7^ since the corresponding values differ in 
the thermal state [3] even for these spherical entangling surfaces. Each individually evolves 
monotonically to their final thermal values as a function of t_4. 

(ii). Strip-like region: Our second example is a strip like region A C M'^"^ in Poincare AdS; 
we take ^ to be a segment of the real line in one of the directions (call it x) and translationally 
invariant in all other spatial directions. The problem of computing the causal wedge in this 
case still continues to be an effective three dimensional problem. In this case it is known that 
Xa ^A even in the vacuum AdSf^+i geometry [3, 7], so this makes for a good example to 
illustrate the general features we argued for in §2.1. We anticipate as described there that xa 
will evolve teleologically and numerical checks show that it indeed does so. Note that while it 
is still true that in the analog of regime 2 (t^ < 0) the surface continues to lie on 9J+[g^], 
the fact that 2^ lies outside in AdS^+i implies that the generators of dJ'^[q'^] have positive 
expansion towards the boundary. With the introduction of the shell, S_4 bends down and toward 
the boundary, thus moving in the direction of the expanding generators and thereby ends up 
having greater area consistent with the general expectations outlined in §2.1. 

5 Discussion 

We have explored properties of bulk causal constructs of [3] which are naturally associated with 
a specified spatial region A in the boundary field theory. In particular, we studied how the causal 
wedge 4^, causal information surface S^, and the causal holographic information xa behave in 
time- dependent bulk geometries, in order to glean further hints for what these constructs might 
correspond to in the field theory. 

While we do not yet have the answer to this important question, and therefore there might be 
no apparent motivation for field theorists to study these constructs, it is useful to draw a lesson 
from the story of entanglement entropy: the work on understanding properties of co-dimension 
two extremal surfaces could likewise have been largely ignored were it not for the connection 




(4.5) 
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with an important field theory quantity, the entanglement entropy; yet most of our insight into 
the behavior of entanglement entropy and related quantities derives from understanding the 
behavior of the bulk surfaces. Certain important properties, such as strong subadditivity, which 
are difficult to prove directly on the field theory side, become manifest in the bulk description. But 
suppose that we did not have this connection between entanglement entropy and bulk extremal 
surfaces yet. Nevertheless, results about extremal surfaces would secretly contain important 
insights, waiting to be realized, about the field theory. While extremal surfaces enjoy the 
status of having an associated field theory quantity Sj^ already identified, our causal constructs 
^A, "^Ay and xa fall into the latter category: We do not yet know what field theory quantity 
they correspond to. We study their properties to gather hints about what such a dual quantity 
might be, but we do not offer a definitive answer. Nevertheless, the results we uncover may bear 
more fruit later when the dual of x and S are finally identified. 

The path towards elucidating the nature of x and H which we chose to follow in the present 
work focuses on a specific class of spacetimes, namely the thin-shell Vaidya-AdS geometries, 
which describe maximally rapid collapse from pure AdS to a black hole. This choice not only 
made our calculations tractable, but offered explicit results in case of physical interest which is 
in some sense furthest removed from equilibrium. Before reviewing these results, let us remark 
on one potential drawback to this approach. In particular, the large amount of symmetry which 
rendered the calculation tractable simultaneously renders such cases somewhat special, so that 
greater caution is needed in drawing general lessons. We have however exercised this caution and 
explicitly identified where and why the calculation simplifies. Moreover, these cases also enjoy 
an important physical significance, as already noted in [3]. 

In the 3-dimensional case, where the bulk geometry is pure AdSs before/inside the shell and 
BTZ after/outside, the causal information surface S_4 coincides with the extremal surface 
only in Regimes 1 and 4 (identified in Table 1) when t_4 < — or > (y9_4, respectively, which 
implies that xa = Sa in these regimes. Fig. 3 explicitly demonstrates that and differ in 
Regimes 2 and 3, i.e. for — v?^ < ^a < '^A- On general grounds, we might then have expected that 
Xa > in these regimes. Nevertheless, we have seen that in fact xa = in Regime 2 as well, 
namely when < 0. This is because the corresponding causal information surface lies on the 
same light cone as CB_4 whose null generators have zero expansion. The important implication of 
this result is that while we might expect on causal grounds that xa behaves quasi-teleologically 
in t_4, in the present case it is completely causal: by measuring xa^ one cannot determine the 
presence of the shell until the shell has been injected on the boundary. 

On the other hand, in Regime 3, where > 0, H_4 is necessarily kinked by the shell. (The 
extremal surface is likewise kinked by the shell, but it no longer lies on the same light cone 
as H_4.) Here we indeed confirm that xa > (cf. Fig. 8), and in the process discover another 
surprise: It is easier to find x than to find S\ Ordinarily one would have expected that due 
to the temporally non-local nature of the causal wedge, finding x is more complicated, since 
it requires us to know the spacetime sufficiently far into the future and past of t^, whereas 
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once is determined, changing the metric to its future or past does not affect it. (In fact, 
direct numerical computation of x is indeed more involved than that of S.) However, in this 
case, lies piecewise either on (future) light cone in AdS or (past) light cone in BTZ, and we 
can therefore evaluate its length by computing the length of different curves on the same light 
cones, connected by flows along non-expanding generators. In particular, the computation of x 
then reduces to finding affine parameter along spacelike = geodesies in pure AdS or BTZ, 
plus finding the intersection of S with the shell, which yields the simple result (3.19). On the 
other hand, computation of S requires finding a geodesic which refracts through the shell, and 
is moreover posed boundary-value problem. 

The analytic simplification in evaluating x allows us to find its scaling behaviour in the 
regions near its initial (vacuum) and final (thermal) values. We find that both x and 5* start 
growing quadratically at small t^, with cubic correction for x and quartic correction for S. On 
the other hand, near the saturation point = y^^, both x and 5* exhibit a faster 3/2 power- 
law behaviour, with greater subleading correction to S than to x- It would be interesting to 
understand the significance of these exponents from a field theory perspective (especially for S). 
In our 3-dimensional setting, we can in fact compare S and x directly, since both have the same 
divergence structure. Evaluating x — 5* (presented in Fig. 8), we find that not only the divergent 
parts, but in fact both the leading finite piece and the first subleading pieces cancel, so that 
X — 5* grows only as t\ when tjs, — ?■ 0+ and as {ifjs, — t^)^ when tjs, — ?■ It is quite remarkable 
that despite the geometric differences in the construction of S_4 and S^, their lengths xa and 
S'_4 agree to such high order. We believe that this is related to the point discussed below of why 
these cases are so special. 

The fact that the deviation from the vacuum value near = is slower than the deviation 
from the thermal value near tjs, = ifjs, indicates where the effect of the shell is greatest. Near 
t_4 = 0, both S_4 and €^ cross the shell near the boundary, where the shell is weak. As t_4 
increases, these surfaces cross the shell deeper, where it gains more strength. Near = 
only the tips of and (i.e. their parts at small if < <^ 7t/2) feel the shell. Although 
this constitutes a tiny region of the full surface which is appreciably affected by the shell, the 
intersection is radially deeper, and the effect of the shell correspondingly stronger. The latter 
effect is the more important one, causing x and S to deviate from their static values more quickly 
near the thermal end than near the vacuum end. 

Having enjoyed the simplifications specified above in three dimensional setting of §3, ulti- 
mately afforded by the fact that and CB^ coincide both in AdSs as well as in BTZ, we briefly 
considered higher dimensional situations in §4. There one of the simplifications disappears but 
we can still consider special cases where another simplification prevails. In particular, for spher- 
ical regions A, the bulk causal wedge 4^ = J^[q^] H J^lq"^] is generated simply, and S_4 still 
coincides with €^ in pure AdS^+i for any d. This implies that in Regime 2, H_4 is still on the 

Although the strength of the shell (e.g. as measured by the amount of refraction of geodesies which cross 
it) blueshifts into the bulk, this might be offset by the portion of the curve which traverses the other side of the 
shell. 
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same light cone in AdS as (along which the null generators always have zero expansion as 
shown in [7]), so that xa = Sj[ for < 0, just as in the 3-dimensional case. However, it is 
no longer true that the extremal surface coincides with the causal information surface S_4 in 
the higher-dimensional black hole, Schwarzschild-AdS(i+i for d > 2. This means not only that 
Xa 7^ in Regime 4, but also that we can no longer find xa as simply in Regime 3. In partic- 
ular, while it is still true that x is piecewise either on a (future) light cone from the boundary in 
AdS or a (past) light cone from the boundary in Schwarzschild-AdS, we can no longer evaluate 
the area of the latter by the same trick of simply comparing with the corresponding part of the 
extremal surface. 

To take another step towards the generic situation, in last part of §4 we drop the other 
simplification as well by considering non-spherical entangling regions dA, while nevertheless 
retaining tractability of the computation. Specifically, in the case of the strip in planar Vaidya- 
AdSrf+i for d > 2, we find that H and (£ do not coincide even for pure AdS^+i, as already noted 
in [3]. In this case we confirm that xa differs from its AdS value already in Regime 2 - in 
other words, here xa does behave quasi-teleologically. As pointed out in §2.1, this is a very mild 
form of teleology: xa knows about the shell only short time before the shell (on the timescale 
of light-crossing transversely across the strip). Nevertheless, this is an important data point to 
keep in mind when searching for plausible field theory duals to this construct. 

As pointed out in [3], one lesson of recent findings is that xa and Sa coincide whenever 
the degrees of freedom in A are maximally entangled with those in its complement A"^, namely 
when the region ^ is a spherical ball in planar or global AdS in all dimensions (or an arc of the 
boundary circle in the BTZ spacetime). The fact that the coincidence extends slightly beyond 
these stationary cases is suggestive. We believe that as a result of this maximal entanglement 
in the AdS geometry for the said regions, it is impossible for xa to grow from its vacuum value 
until the disturbance has come to play on the boundary. It is for this reason in the field theory 
that despite the deformation of the causal wedge in a quasi-teleological fashion, xa nevertheless 
evolves causally. 

So far, we have restricted our attention to a specific class of configurations, namely the thin- 
shell Vaidya-AdS spacetimes, serving a convenient toy model of quantum quench in the field 
theory. A complementary approach to elucidating the nature of x and S is to maximally relax 
the assumptions about the bulk spacetime, and consider general global properties that these 
constructs must satisfy. Since this approach is rather more formal, we have chosen to present 
the results separately in a companion paper [45] . As a preview, here we simply note some of the 
key features. 

It is sometimes useful to consider disconnected regions, A = AiU A2 with ^1 fl ^2 = 0- If 
Ai and A2 are taken to lie at the same time in the boundary field theory, then it follows trivially 
that the causal wedges for the two parts are disjoint, and the causal holographic information is 
simply the sum of the two parts. More generally (for intersecting regions) it is easy to establish 
subadditivity, though as demonstrated by explicit examples in [3], strong subadditivity is not 
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universally satisfied. In higher dimensions we can also consider a single but non-simply-connected 
region A. In such a case, 4^ is likewise non-simply-connected, and consists of multiple 
components. More curiously, even for simply-connected regions A, 4^ itself need not be simply- 
connected, and S_4 may consist of arbitrarily many disconnected components. We explain and 
explicitly demonstrate this in [45]. 

While the above remarks might lead the reader to expect that "anything goes" and that 
causal wedges have properties which are hard to characterize globally, there are some features 
which hold in full generality. One such important feature, already alluded to above, is that causal 
wedges can never penetrate the event horizon of a black hole. This follows simply from causality: 
causal curves from within the black hole can never reach the boundary and therefore the interior 
of a black hole cannot be contained in the causal wedge of any boundary region. Nevertheless, 
♦_4 can reach close (or up to) the horizon for suitable regions A. 

In this context, there is an important difference between global and planar asymptotically 
AdS geometries: For the field theory living on Minkowski background, there is no upper bound 
on the size of A, and S_4 can penetrate arbitrarily close to the horizon for arbitrarily large regions 
A. In fact, in this regime it is easy to see that the finite piece of xa scales extensively with the 
volume of A. On the other hand, for the field theory living on Einstein Static Universe, the 
region A can at best wrap the sphere. That means that the extent of the domain of dependence 
is either bounded by ip_A_ < tt or fills up the entire boundary spacetime. In the former case, 
only reaches a finite distance from the black hole. On the other hand, in the latter case, (for 
both ESU and Minkowski boundary geometries) when A covers the entire Cauchy slice of the 
spacetime, the boundary of the causal wedge by definition coincides with the event horizon. In 
this case xa is precisely the black hole entropy. 

Finally, let us contrast this feature of the causal information surface with the extremal 
surface While the black hole deforms outward as compared to the AdS case (with same 
ipj\), and in static bulk black hole spacetime extremal surfaces anchored on the boundary must 
lie outside the horizon (both of these features were recently demonstrated by e.g. [46]), in a time- 
dependent situation the extremal surface can actually penetrate the event horizon. This was 
argued already in [46, 48], and seen explicitly in [17] for the planar Vaidya-AdS case, but is also 
manifest in the right panel of Fig. 3. These issues are examined further in the forthcoming work 
[47]. Hence while causal holographic information is not cognisant of the causally disconnected 
region inside a bulk black hole, the entanglement entropy does sample at least a bit of the 
spacetime inside. 
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A Computational details for in Vaidya-AdSs 

In this appendix wc present some details about the computation of quantities which are relevant 
for the explicit determination of S_4 in thin shell geometries. 

Let us begin by collecting expressions of the null geodesies occurring in the determination 
of for —'•PA<'tA<VA- III order to shorten the formulae, we find it convenient to introduce 
V^ir) =T^- (L/E^fUr); explicitly 

V\r) = {1- - e , V°{r) = (1 - fy + frl . (A.l) 

In the following, these quantities will occur with a further subindex which is either y> or -x, 
indicating whether the corresponding geodesic is respectively outgoing or ingoing. 

• An outgoing geodesic lying entirely outside the shell and arriving at reads 

vJf, 7+) = — log ; , f 7 r + (A. 2) 

• Likewise, an outgoing geodesic lying inside the shell and connecting to a point of the 
shell with coordinates Ps = (^s = 0, r^, (ps) is given by 

v^(r, = arctan yjv),(r) + arctan(r) — arctan yjv),{rs) — arctan(r5) (A.4) 

ipy,{r,k+) ^ arctan (fe^^V),(r)) - arctan (^^V),(rs)^ + (A.5) 
The choice of the integration constants guarantee that v ^{ts, k+) — and (p ^{r, k+) — (ps- 



An ingoing geodesic outside the shell which connects a point of the shell to a point of 5_4 
reads 

^^r^J-) = 7^ log /T7^ zrr^ - 7^ log /YTw^ zr-r^ (A-6) 




^K°(r) -Th r + ThJ 2rh \^/V^(rsj - Th rs + rn 
= — log ^7=^^ : /tt^^ • + "Ps (A.7) 
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The integration constants are obtained by imposing v^X^syj-) = and V9^(rs,j_) = ip, 



• Finally, an ingoing geodesic which starts at is given by 

f^^(r, = fv — arctan ^^K((r) j + arctan(r) (A.8) 
ip^{r,k_) = - arctan + |sign(fc_) (A.9) 

A.l Solving the refraction conditions 

Having explicit solutions for the geodesies we turn to the refraction conditions. We discuss an 
alternative way to find the solutions of the refraction conditions across the shell which is slightly 
different with respect to the one described in §2.2 (and serves to provide a complementary 
viewpoint). The equations to solve are 

1 1 _ l + rl viir,) <(r,) ^^^^^^ 



Given a geodesic crossing the shell at the point ps from one side, these equations tell us which 
is the corresponding geodesic on the other side of the shell. We find it convenient first to find 
the solution from the second equation of (A. 10) and, subsequently, employ the first one as a 
consistency check. 

In order to deal with the second equation of (A. 10) notice that 

The second equation of (A. 10) can be written as Ci{rs,k) = Co{rs,j). From this equation we 
can extract either k{rs,j) or j{rs, k) finding 

^ ^ 2rsCo{r„j) ^ ^jrsfojrs) 

C!{rs,j)Mrs) + l r,[/,(r,) + /„(r,)]+r/[/,(r,)-/„(r,)]0^ 
. ^ 2rsQ{r,„k) ^ 2krji{rs) 



The second expression in (A. 12) and (A. 13) is obtained multiplying the first one respectively 
^ ^ _ rs-v^yv"(rs) ^ ^ rs-nyJvKrs) _ ^^^.^^ ^^^^ .^^^2) and (A.13) cau be interchanged by 

exchanging the inside and outside quantities, as expected. 



A.2 Regime 2: -y?^ < < 

In this regime the whole refraction curve belongs to 9+ (4^) and therefore lies entirely inside 
the shell. 
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A. 2.1 Refraction curve on 9+ (4^) 



The outgoing geodesic arriving at crosses the shell at the point ps, whose radial coordinate Vg 
is defined by v /.(rg, j+) = 0. From (A. 2), this equation reads 



e-'"'^^^ . (A.14) 



For the radial case j+ = 0, the l.h.s. of (A.14) simplifies to a square, giving = rh coth.{rhV/s^/2). 
In the generic case, writing (A.14) in a form where ^/V°Jj^ is isolated on one side of the equation 
and then squaring it, we obtain an algebraic equation of fourth order in which admits = ±1 
as solutions. Thus, we are left with the following second order equation 

9 , / N 1 — ?i coth^fr/itiA) r, , 
tI - 2rh co\\y{rnV^) + _ = , (A. 15) 

whose largest root is the first formula in (3.8). The angular coordinate = V'/'(^S5J+) can be 
found from (A. 3). By observing that V^ir^ reduces to a perfect square, one obtains the second 
formula of (3.8). 

The outgoing geodesic of c}+(#yi) which refracts at the point of the shell is made by the part 
inside and the part outside the shell, which are characterized by /c+ and j+ respectively. These 
coefficients are related one to each other by (A. 12) and (A. 13) with 77 = +1 (see e.g. (3.9)). 

A. 2. 2 inside the shell 

In the regime we are considering S_4 is entirely inside the shell. The curve is the solution of 
(3.10). In the first equation, the term arctan(r^) cancels because it occurs both in (A. 4) and 
(A. 8). In order to deal with (3.10), first we bring the terms dependent on r* on one side of the 
equations, then we take the tan of them. Employing the addition formula tan(a + 6) = i^^^^^^^ih 
and the property tan(x ± 7r/2) = — cot(x), we find 



^^^^/Et ^t!! ^ = f^v + arctan .M(r,) + arctan(r,)) = (A.16) 

1 - vv^>(n)v:((n) ^ V y 

V^5pit±^M^ = _ eot (arctan Ll^)lk\ - ^^ ^ (A.17) 

In the special case of the radial geodesies k+ = k- = 0, the equation (A.16) simplifies to 

2r= 

tan [^Vy + 2arctan(rs)) . (A. 18) 



1-ri 



Then, using the duplication formula for tan in this equation, we find the solution (3.11). 

For non radial geodesies, we can obtain simpler expressions for (A.16) and (A.17) multiplying 

their l.h.s. 's by 1 = "/^-^^ ^^"^^ and 1 = ^^-^^ "^'^ ^r-^^^h-^ , — ■, respectively. Besides the 
radial case, this step is not allowed for kj^ = k-, but this is never realized in this regime. This 
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trick leads to the important simplification of a factor (r^)^ + 1, allowing us to write (A. 16) and 
(A. 17) as a linear system in terms of \fV~y<\ 



and a/K 



K^^r^), whose solution reads 



(1 

+ fc_)(A;+T; 



{l-kl 

Taking the square of these two equations, we can find 



V (p 



(A. 19) 
(A.20) 



J\2 



and a consistency condition 



1 



k' 



(k^ + k_){Ti-k^T^ 



'i-ki 



ki + 



{k^ + k.){Tl-k.Tl 



kl 



T T 



'~ri'~ri 



1^ 



kl + 



ik+ + k.){Ti-k+T^ 



kl 



'T'i'T'i 

V if 



(A.21) 



(A.22) 



In order to understand these relations, we recall that (r^, ips) and fc+ depend on {j+,V/^) through 
(3.8) and (3.9). Then, from (A.16) and (A.17) we write = T^(j+, t;^, ^v) and = T^{j+,V/,). 
Plugging these dependences into (A.22), and inverting (numerically) this relation we obtain A;_ = 
k-{j^, fA, fv), telling us which geodesic of ^-(♦^) intersects the geodesic of 9+(#^) characterized 
by j+. Substituting this result into (A.21) we finally find = r* (j+, f^, fv)- The angular 
coordinate (p]^ = ip^{r]^, /c+) = (y9^(j+, f a, fv) is then obtained through (A. 5). 
As a check of (A.21) and (A.22), notice that in this regime of tjs, and Vj+ we have < 0, namely 
that S_4 is entirely inside the shell. 



A. 3 Regime 3: Q < tj, < ipj. 

In this regime, the central part of the refraction curve belongs to 9+(4_4). The corresponding 
part of H_4 is inside the shell and it can be found as explained in §A.2.2. Thus, the results found 
above must be applied only for < < 1 or, equivalently for |A;_| < A;* < 1, where j'^ (or 
k*J) characterizes the critical geodesies. 

For < < 1 (or A;l < |A;_| < 1 equivalently), the refraction curve belongs to (9_(^_4) and 
therefore the corresponding parts of S_4 (they are symmetric w.r.t. to the radial geodesies) are 
outside the shell. For these geodesies, the results showed in this subsection must be used in order 
to find outside the shell. 

A. 3.1 Refraction curve on 9 (4^) and critical geodesies 

The ingoing geodesic starting from at the boundary crosses the shell at ps first. The radial 
coordinate of this point is such that f^(rs, k^) =0, i.e. from (A. 8) 

arctan (^^Jv^tjjs)^ ~ arctan(rs) = t>v . (A. 23) 

First, we take the tan of this equation employing also the subtraction formula for the tan. Then, 
multiplying the l.h.s. of the resulting equation by 1 = ^, — -, we find expression where 
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\/VHrs) occurs only in the denominator. Isolating \/VHjs) and then taking the square, we 
obtain the following second order equation for 

rl - 2 cot(t;v) -_ ^ = , (A.24) 

whose largest solution is given in the first equation of (3.12). The angular coordinate ips = 
V?K^(rs,/c_) given in the second equation of (3.12) can be computed through (A. 9), by further 
noticing that K^(rs) becomes a perfect square. 

The ingoing geodesic starting from g^, characterized by /c_, becomes an ingoing geodesic char- 
acterized by j_ and propagating inside the shell. The function j_{kJ) is given in (3.13) and it is 
obtained from (A. 13) with rj = —1. 

In this regime of t^, two symmetric critical geodesies on ^-(♦^) occur. They meet both the shell 
and d^{^X) the same point. This implies that the corresponding critical value fcl satisfies 
(A.24) with given in (A. 21) instead of r^, namely 

. , x9 . X 1 + (k* )^cot2(i;v) , , . , 

«)2-2cot(t;v)< i_(^k*y ° ^ ^^■'^^^ 

As discussed in the end of §A.2.2, from (A. 21) and (A. 22) we find kZ_ = k*_{j\,t_A,ip_A) and 
= ''^x(i+5 V^^)- Substituting these results into (A. 25), it becomes an equation for i\ that 
we can (numerically) invert, getting j*^ = j!^(t^, fA)y where < t_4 < (/j^. 

A. 3. 2 outside the shell 

The curve S_4 outside the shell is given by (r°, ip") satisfying 

v^{r:,j.) = vAr:,u) = < , V^{r:,J-) = ^ACu) = (A.26) 

where the geodesies of 9+(#^) are given by (A. 2) and (A. 3), while the ones belonging to (9_(#^) 
are described by (A. 6) and (A. 7). Notice that in the first equation of (A.26) the term ^ log 
simplifies. Writing (A.26) as equations involving the arguments of the log's, they become 



(A.28) 



Muhiplying the l.h.s. of (A.27) by 1 = ^^^^ and its r.h.s. by 1 = ^^E^, the 
equation simplifies. A similar simplification occurs in (A.28) when we multiply its l.h.s. by 
1 = "^,5! ~ — - and its r.h.s. by 1 = ""l resulting equations, the dependence 

of r° can be isolated on one side, which becomes the square of a simple rational function in terms 
of ^JV°{r'^■ Taking the square root of the two equations in such form, they can be written 
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respectively as follows 



v/K°(r°)+rfe ^ 
V^l/^(r°)+j_rfe 



^(l-J^)(r.-rO(v/T^^ + r,) 
(1 - j2 )(r, + r^) (v/V^(^ - ry) 



rpo 



(A.29) 
(A.30) 



It is now clear that these equations can be written as a linear system in y^V°{r°) and \/V^(r°), 
which can be easily inverted, giving 



t^ -j.t: + {i-u)t:t^ 

V (p 



(A.31) 
(A.32) 



Squaring both these equations, we obtain (r°) and a consistency condition 



•A 



i^j- + T;-]+T' 



'To "To 

V ip 



7^0 _ 7^0 



t"-j.t: + {i-u)t:t^ 



'~Po 7^0 

V ip 



(A.33) 

(A.34) 



These equations allow to find the part of S_4 outside the shell, as well as (A.33) and (A.34) lead 
to the part of S_4 inside the shell. Indeed, first we observe that (r^, v?s) and j_ depend on (fc_, fv) 
through (3.12) and (3.13). Then, the definitions (A.29) and (A.30) tell us T° = T°{k_,Vy,v^) 
and T° = T°(A;_,iv)- By inserting all these functions into the consistency condition (A.34), it 
becomes an equation which gives us (through numerical inversion) j+ = j+(A;„, fv, f a), namely 
the geodesic of 9+(#yi) intersecting the geodesic of (9_(#^) characterized by Given this 
result, (A.33) allows to obtain r° = r°(fc_, fv, f a)- The angular coordinate ^9° = y9,^(r°,j_) = 
(/9°(/c_, fv, f a) is then obtained through (A. 7). 

We can check (A.33) and (A.34) by verifying that in this regime of tj^ and for |A;_| > /cl we have 
f ° > 0, namely that is outside the shell. 



B Computational details for Sj^^ 

In this appendix we give some details for the analytical computation of S*^ in three dimensions. 
We generalize the discussion of [20], allowing for the geometry inside the shell to be either 
Poincare AdS or global AdS or another BTZ spacetime (corresponding to heating up a preexisting 
thermal state). While many of the results follow from our discussion in the main text and 
Appendix A modified appropriately to spacelike geodesies, it is useful to work these out explicitly 
to obtain compact expressions for the entanglement entropy. 
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B.l Spacelike geodesies and refraction conditions 



While for the determination of S_4 we needed the null geodesies, for the holographic entangle- 
ment entropy 5*^ we have to compute the spacelike geodesies anchored on the boundary at the 
endpoints of the interval A. 

We can be sufficiently general and consider a metric like (2.1) with d = 2 and f{r,v) = 
fi{r) + 9(f)[/o(r) — /j(r)], where we always consider a BTZ geometry outside the shell (i.e. 
/o(^) = — r^^), while inside the shell we choose either Poincare AdS (/i(r) = r^) or global 
AdS (/i(r) = + 1) or another BTZ (/i(r) = — r^^J with Vh^i < rh,o in order to satisfy the 
null energy condition. 

As discussed in §2.2, the spacelike geodesies (k = 1 in (2.9)) are characterized by the pair 
{E,L). The ones we are interested in are made by three branches: two disconnected and sym- 
metric ones outside the shell and one inside the shell. The coordinates of the meeting points 
are {vg = 0,rs,±ips). The symmetry of the problem allows us to consider only a half of the 
geodesic, the outgoing one, going from (f(£,r£,0) to the point (t^, 00,^9^) of the boundary. Be- 
ing (fi continuos across the shell, Li = Lo, while the jump of E is given by (2.16), where k = 1 
occurs through i) at the shell (see (2.14)). Since we are considering equal-time endpoints (for 
endpoints at different times on the boundary see [44]), inside the shell we have {Ei, Li) = (0, r^). 
Then, from (2.14) we find that v at the shell for spacelike and outgoing geodesies {1] = 1) reads 
"^s = [{r'^ ~ ^'i) / {i^sfii''^s))Y^'^ and, through (2.16), this leads to Eq. Thus, the pair {Eo, Lo) reads 

[fo {rs ) - Mrs) ] ^rl - r\ 
2rsA/^ 



Eo = J^vsJ^y^s 'e ^ Lo = Li = r^. (B.i; 



The equations to solve in order to find the spacelike geodesies are 

t' = , , ^' = , = , (B.2) 

faVEl + fa{l-Ljr^) r^^El + Ul-Ljr^) 

where a is either i or o, and fi{r) is one of the three choices described right above. 
Outside the shell we have a BTZ geometry and the solutions of (B.2) are 



2 



U{r) = — log \+C,^ t{r) + Ct (B.3) 




{rh + Eo)rh + V^/Djr) 



where we introduced Do{r) = E'^r'^ + (r^ — T^)(r^ — L^). It is important to remark that the 
integration constants Ct and C^ provide the boundary data, namely 

Ct = t^, C^ = V^. (B.5) 

These parameters are also related to (r^, (ps) of the point at the shell by the conditions Vo{rs) = 
and (pi{rs) = ipo{rs) = (ps, which give respectively 

tA = - i{rs) - J-log( ^' ) , 99^ = - <^(r^) . (B.6) 



2r/i \rs + rh 
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Let us focus on the first equation of (B.6) and try to extract r^. First, notice that, from (B.l), 
we find that -Do(^s) = {r^ - rl)[fo{rs) + fi{rs)f /[4:fi{rs)]. This result gives 



1 



2(^2 - rl)r, + [rs\fo{rs) + /.(r.)| - r,(/,(r,) - /,(r,))] 



t(r.) = ^\og h==^ ^ ■ (B-7) 



2rh 



2(^2 _ ^2)^^ + [r,|/„(r,) + /,(r,)| - r,(/,(r,) + /,(r,))] 



Following [20], we find it convenient to introduce a new parameter 6 which mixes and r^. Its 
definition depends on the geometry inside the shell 



cose ^ J = ^ - + + 1) global AdS (B.8) 

Notice that the cases of Poincare AdS and global AdS are obtained from the expression for BTZ 
by substituting respectively ^ = — 1 and ^ = 0. This feature persists also in the expressions 
hereafter, unless otherwise specified. From (B.8), notice that can be written in terms of 
and 6^ as r| = — (rf — •) cos^ = rl sin^ 6 + • cos^ 6. 

Plugging the definition (B.8) into (B.7) with the assumption fo{i^s) + fii^s) > 0, the first equation 
of (B.6) leads us to write as a function of tyi and 6. After some algebra, one finds the following 
algebraic equation of the second order in 

—) ' - coth(r,,„t^) ^ + f 1 - = (B.9) 

rh,oJ rh,o l + cos^ y rl^J 

whose largest root reads 



a/I — r|/r 2 Poincare AdS 



rho 



coth(r,,„t^) + ^/coth^(r,,,t^) - ^""^^^^^^ jf/"^"'"^ I ■ (^-10) 



Notice that is contained also in the r.h.s. of (B.IO) through cos^, and this means that we 
have not fully inverted the first equation of (B.6). Nevertheless, the formula (B.IO) gives in 
terms of and 6. 

From the second equation in (B.6) we can write (fji^ in a more compact form. From (B.4) with 
rj = +1 and Do{rs) given above, we obtain 



1 / 2(1 + cos 6')r^ - 2r/jo./r2sin^6' + -cos^e^ - {jl ^ + rl ^) cos 9 

^W = — log| I ' = ■ ' I , (B.ll) 

2(1 + cos 6')r2 + 2rh,o\ r1 sin^ 6 + • cos^ — {r\^ + r\ ■) cosO 



2r 



h,o 



where the term with the square root is simply extracted from (B.8). 

In order to find (fs = fi{fs), we need the solution ipi{r) of the second equation in (B.2) for the 
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different ctioices of /j. Tliey read 



r,f r 




1 — r^)r^ — 2r|\ tt 
(l + rl)r^ J ^ 4 



r2 + 



r;,,j + ^(r^ -r^ .)(r2 -r|) 



Poincare AdS 
global AdS 

BTZ 



(B.12) 



wliere tlie integration constants liave been fixed by imposing tliat (pi{r^) = 0. Now we can write 
fii'Ts) = in terms of and 9 as follows 

f cote 



arctan 



COS0 



a/ sin^ 6 — cos^ 6 




rl sin^ ^ + "^hi ^os^ ^ + '^h,i cos 6* 



Poincare AdS 
global AdS 

BTZ 



(B.13) 



sin^ ^ + j cos^ — r^^i cos 6* 



For 6* G (0,7r/2), the first expression is the limit r^^i — )■ of the third one, as expected. As for 
the second one, it can obtained by setting rh,i = i (the imaginary unit) in the third one and then 
using that ^ log(2;/z) = arg(z) for any complex number z. 
At this point we remark that 



(B.14) 



The first equality comes from the second formula in (B.6), (B.ll) and the expression of cps in 
(B.13) corresponding to the geometry inside the shell, while the last step is found by further 
substituting = rs{tjs,,9) given in (B.IO). Thus, by (numerically) inverting (B.14), one can find 

= 0{tA,^A)- 

Behaviour near = and tj, = ^px- It is instructive to perform such inversion through 
series expansions around = and tj, = ^j^. In both these cases we use that 

tan V5s(tyi, 6') — tan <^(/f:^, ff) 



tan ipj^ 



(B.15) 



1 + tan i-PsitA-, tan V9(t^, Q) 
where ^s{^a-, ^) is the expression in (B.13) corresponding to BTZ and <^(t^, &) = '^{r^ is obtained 
by plugging (B.IO) into (B.ll). We stress that all the following results provide the corresponding 
quantities also for Poincare AdS and global AdS inside the shell. They are easily obtained as 
the special cases of r/^ j = % (the imaginary unit) and r/^ j = respectively. 

About = 0, using that 6'(0, ^ = 0, we can write 

+00 +00 

0{tA,VA) = J2^''^^'^^^A ^ tan^^ = ^$„((^^)t:^ (B.16) 



k=l 



n=0 
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where for tany?^ (B.15) is applied and therefore the coefficients $„'s depend on Lpj^ through the 
Gfc's. Since the Lh.s. of the second equation in (B.16) is independent of t^, it can be solved 
order by order, namely tan(y9_4 = $o(v^y4) and = $„((y9_4) for n ^ 1. From these equations we 
find respectively 

61 = —-^^^^ r, 62 = 0, 63 = ^ ^ ^, 94 = 0, ... 

smh(rft.j ^pX) 12 

(B.17) 

which give the expansion of around tj^ = 0. 

As for the expansion of 0{t_A,ipj\) around t_A = notice that 0{ipjs,,({>jC) = 7^/2. Now we 
need to introduce half-integer powers, namely (recall that < < 9?^) 
+00 +00 

^(U, ^a) = ^-Y^ Qk/2{^A) {^A - tA?'^ =^ iailipA = ^p/2{^a) {^A - tAf'^ ■ 

k=l p=0 

(B.18) 

As above, we can solve the second equation order by order. The first order gives no information 
since from the expansion of the r.h.s. of (B.15) with the first formula of (B.18) plugged in we 
find $o(<y'^) = tanyj^. For p > we impose $p/2(v'^) = 0, finding that 

2 ^2 



61/2 = \ 2rh,oCot\i{rh,oVA) , ©i = — — tanh(r/j,o v^^) , 

V 3 rh^o 

r, -2 (5 " 4 rl- „ + 2 J + 3 „ csch2(r;,,„ ^j,) + 4 (r^ „ - sech2(r^,„ ^j,) 

"3/2 — . 

18 rl o V2 Th^o coth(r,,,o ^pX) 

(B.19) 

and so on, which provide the expansion of Oitji^.pX) around = pji,. 

In order to find e.g. = r^itji^, pX)i we first write r\ = rl — {rl — J cos 6 from (B.8), then 
we plug in rs{t_A,6) given in (B.IO) and finally we employ 6 = 0{tji^,{pjC) obtained by inverting 
(B.14), as explained above. A plot of = r(c(t_4, pX) as a function of for a fixed is shown 
in Fig. 6. Through the expansions of 6{tj[, pX) found above, we can write the first terms of the 
expansion of rg = r^itj^, ipj) when — )■ and — )■ V9_4 for a BTZ spacetime inside the shell. 
They read respectively 

u ^ fi.^ , , irjo - rl^yji coihjrh^i p^ ,4 , .0 o^^ 

r^tA^ '^Aj = rh,i coth(r;,,i pX H „ . , 2/ ^ t_^ + ... (B.20) 

8 smh {rh^i pX 

tSa^^a) = rh,oCoth{rh,o^A) /o ^u/ ~ ^ + • • • • 

^2r/,,o coth(r/,,o V'^) Sr^^^coth {rh,oVA) 

(B.21) 

We can also find in terms of (fu,ru) at the turning point. Indeed, inside the shell we have 

Vi{r) = U{r) + p,{r) + Bi (B.22) 

where Bi is a constant and Pi{r) is the solution of p' = l//i(r) given by Pi{r) = log 

for BTZ, Pi{r) = arctanr for global AdS and Pi{r) = — 1/r for Poincare AdS. Since Ei = 0, we 
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have that tj(r) = for every r inside. In particular, considering the point r = r^, we obtain 
Bi = v<f — Pi{r^). Then, by imposing Vi{rs) = 0, we find that Pi{rs) = Pi^r^) — 1>£ and therefore 
we can write in terms of (fs,^-^). 



B.2 Holographic entanglement entropy 

According to the prescription of [5-7], the holographic entanglement entropy S*^ is given by the 
length of the spacelike geodesic studied in the previous subsection. The functional giving the 
length of a piece of curve which is either entirely inside {a = i) or entirely outside the shell 
(a = o) reads 



= / \ -Ut'Y + ^ + r^ip'f dr . (B.23) 



Integrating (B.23) with (B.2) plugged in, a primitive is 

£„(r) = ^ log (e'^ - + 2r2 - + 2^ E^r'^ + {r'^ - rl){r'^ - L"^)^ . (B.24) 

We recall that if a = i, then r\ = rf^^ in this expression and, in particular, r^^ = — 1 for global 
AdS and r^ - = for Poincare AdS inside the shell. In the latter case the argument of the 
logarithm becomes a perfect square. 
The holographic entanglement entropy is given by 

AG%^ Sa = 2[A(r,) - Ciir^)] + 2[C^^^ir^) - C^^^ir,)] (B.25) 



where G^^^ is the three dimensional Newton constant and Too — cxd is the UV cutoff in the 

y(3) 



boundary theory (cf.. Footnote 10 for conversion between G)^ and the boundary central charge 

Ceff). 

Outside the shell, by using (B.24) and that £3x2(^^00) = log(2roo) + 0{r^) asymptotically, 
we find 

2[£BTz(roo) - C^Tzirs)] = 21og(2roo) - log - + 2rl - r^ ^ + 2^/D^^ (B.26) 
^ 2 .og(2 - log ( '^"-''-g;-^}^-'^-' - rl + 2ri - + ^iMn) + | 

where Do{r) has been defined after (B.4). In the second step, (B.l) and the expression of Doirg) 
given after (B.6) have been used. Inside the shell, since Ei = 0, the argument of the logarithm 
in (B.24) becomes a perfect square and therefore we find 

2 [A(r.) - A(r,)] = 2 log ( + - rj \ ^^ ^T) 
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Combining together the two contributions (B.26) and (B.27), and discarding the divergent term 
21og(2roo), we find for the regularized holographic entanglement entropy 



rcg 




{rlo - rl,,) 2,02 2 , / 



It is interesting to rewrite this formula by employing (B.8). First, notice that 

VM^s) + \/rl - r\ _ 1 + cos 6* 



2 




(B.29) 



r.2 _ „2 sin Q 

<£ ' h,i 

After some algebra the argument of the logarithm in the second line of (B.28) can be written as 



cose{l -rl -frlj ^ r. 



2{1 + cos 6)rs/rh,o rh^o 



'''"^ ^ (B.30) 



sinh {rh,otA) 



where the second step is obtained by isolating coth{rh,otA) in (B-9) first and then by using 
the identity coth^x — 1 = 1/ sinh^x. Thus, putting together (B.29) and (B.30), one can write 
the regularized holographic entanglement entropy (B.28) in a compact form. Adding also the 
divergent term coming from (B.26), the holographic entanglement entropy reads 

QU \ 1 /2roo sinh(r;,,ot^)\ 

SA{tA,VA) = ^ log ^^77 r B.31 

where we remarked that 9 = 6'(t^,(y9^), as discussed after (B.14). For Poincare AdS inside the 
shell, we recover the result of [19, 20]. Notice that the expression (B.31) is formally the same 
for all the choices of the geometry inside the shell that we considered. The definition (B.8) 
distinguishes among them. 

Behaviour near t_4 = and t_4 = ipy^: It is useful to write the first terms of the expansions 
of Syi{tji^, (/?_4) around = and t_4 = yj^. 

As for the expansion around = 0, first we plug into (B.31) the expansion of 6'(t_4, (y9_4) given 
in (B.16) and then employ the explicit expressions for Qki'^A)^ obtained by solving the second 
equation in (B.16) order by order in (see (B.17)). Thus, we find that 

- SAtA, Va) = 2 log sinh(r,, cp^) + "Jb^qJj^ (B.32) 



(1^2 ^2 / Ci 

' h,o ' h,i I " ' h,i 



+ <o + 3<Jt^ + 0(4 



48 \?,mh^{rh,i(pA) 

When the spacetime inside the shell is global AdS (i.e. when r^^i is the imaginary unit), the 
expansion (B.32) becomes (3.20). Instead, the limit r/^ j — )• of (B.32) provides the corresponding 
result for Poincare AdS inside the shell, namely 

— SAitA, Va) = 2\og{2r^^j,) + '^t\-'j^[ — + rl]t\ + 0{t\) . (B.33) 

Ceff ^ 45 J 
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Discarding the constant term in and taking the hmit of large ^pj, of the remaining part, we get 
— r'lt'^/AS + O(t^). These are the first terms of the expansion for small of the function 
41og(cosh(r/ii(:_4/2)) found in [44] (see also footnote 13). 

A similar procedure provides the expansion of v^^) for tj^ — > ipj^: plugging the expansion 

of 0{tj\^,ipjO defined in (B.18) into (B.31) and then using (B.19), we find 



— SAKtA^^A) = 2 log smh(r;,,o <^^) - ^ , , , , yPA-tA) 

Ceff \rh,o J ^^Jrh,oCoth.[rh^o'^A) 

(rlo-rlif tanh^(r,,,o<^^) 



3rL 



i^A-tA) +••• (B.34) 



where the dots represent subleading contributions. Specializing this result to the case of global 
AdS inside the shell, we find (3.23). Instead, setting rh^i = we obtain the corresponding 
expansion for Poincare AdS inside the shell. 
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